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ABSTRACT 


The  principal  task  of  this  program  was  to  develop  a  better  understanding 
of: 

)  the  mechanisms  of  the  growth  of  passive  films  on  metal  surfaces, 

ii)  the  role  of  space  charge  in  determining  the  properties  of  the  film, 

and 

iii)  passivity  breakdown  phenomena. 

Experimental  data  for  the  ac  impedance  and  passivity  breakdown  properties 
of  passive  films  on  nickel  are  presented  in  Sections  2  and  3.  Also  a 

theoretical  development  of  the  point  defect  model  is  presented  in  Section  4. 

A  comparison  between  the  theoretical  and  experimental  data  for  the  growth  of 
passive  films  on  nickel  and  iron  is  included.  In  addition,  the  role  of  minor 
alloying  elements  (such  as  molybdenum)  on  pitting  and  passivity  breakdown  is 
considered  theoretically,  and  the  model  is  compared  with  experimental  data  for 
18%  Cr  ferritic  stainless  steel. 

The  results  of  this  study  indicate  that  the  point  defect  model  although  at 
present  in  a  crude  state  of  development,  is  a  reasonable  first-approximation 
theory  for  explaining  some  of  the  basic  phenomena  related  to  the  growth  and 
breakdown  of  anodic  passive  films  on  metal  surfaces.  On  the  basis  of  data 
developed  in  this  program  it  is  recommended  that  further  experiments  be 
conducted  to  further  refine  the  model. 
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the  increase  in  the  measured  capacitance  with  potential  probably  arises  from  a 
corresponding  increase  in  the  space  charge  capacitance  of  the  p-type 
0-Ni ( OH ) 2 -  Nevertheless,  a  typical  double  layer  capacitance  also  exhibits  a 
potential -dependence ,  and  part  of  the  observed  increase  in  Cp  may  reflect  this 
factor.  In  addition,  the  relative  contribution  of  the  double  layer 
capacitance  to  the  measured  capacitance  could  well  increase  at  higher 
potentials  where  the  space  charge  capacitance  is  larger. 

Summary 

Impedance  data  for  single  crystal  nickel  in  phosphate  solutions  have  been 
analyzed  in  terms  of  an  equivalent  circuit  model  to  provide  a  parallel 
capacitance,  charge  transfer  resistance,  and  Warburg  coefficient,  o.  The 
value  of  o  was  found  to  be  independent  of  potential  within  experimental  error 
at  pH  values  of  9,  11,  and  12.  At  other  values  of  pH  in  the  range  studied  (7 
to  12),  the  lack  of  a  true  steady  state  for  film  growth  and  a  non-uniform 
composition  of  the  film  may  have  been  the  cause  for  relatively  large  changes 
in  a  with  potential . 

The  presence  of  a  large  charge  transfer  resistance  indicates  that  the 
thickness  of  the  surface  film  is  small  enough  to  allow  both  an  interfacial 
electrochemical  reaction  and  solid  state  mass  transport  processes  to  control 
the  film  growth  kinetics. 

The  capacitance  appears  to  originate  largely  from  the  double  layer  in 
solution,  although  the  space  charge  layer  in  the  B-NilOH)^  film  might  provide 
an  important  contribution  depending  on  the  potential. 


case,  a  rotating  ring  disk  electrode  study  (12)  of  nickel  has  shown  that, 
despite  the  predicted  thermodynamic  stability  of  the  Ni(0H^  ion,  there  is  no 
significant  dissolution  of  the  electrode.  Nevertheless,  significant 
dissolution  at  pH  values  of  7,  and  possibly  8,  may  substantially  reduce  the 
influence  of  solid  state  mass  transport  in  controlling  the  overall  rate  of 
film  growth. 

It  is  difficult  to  determine  the  extent  to  which  the  lack  of  a  true  steady 
state  and  a  non-uniform  composition  affect  the  film  growth  kinetics, 
particularly  solid  state  mass  transport  processes.  Nevertheless,  based  on  the 
drift  in  the  "steady-state"  current,  the  overall  uncertainty  in  and  o  is 
estimated  to  be  ±5*.  It  is  of  significance  to  the  point  defect  model  that, 
when  this  uncertainty  is  taken  into  account,  a  is  found  to  be  approximately 
constant  for  several  pH  values. 

The  presence  of  a  large  charge  transfer  resistance  indicates  that  the  rate 
of  film  growth  is  partly  controlled  by  an  interfacial  electrochemical 
reaction.  At  this  stage,  it  is  possible  to  only  speculate  whether  this 
reaction  occurs  at  the  metal/film  or  film/solution  interface.  Although  it 
appears  from  Table  1,  that  does  not  exhibit  a  well-defined  dependence  on 
pH,  there  is  some  evidence  from  another  study  (8)  that  the  rate-determining 
reaction  may  occur  at  the  film/solution  interface. 

Clearly,  the  film  does  not  become  sufficiently  thick  to  allow  solid  state 
mass  transport  processes  to  dominate  control  of  the  film  growth  rate.  This  is 
consistent  with  previous  work  (9)  which  has  shown  that  the  thickness  of  films 
grown  potentiostatical ly,  or  gal vanostatically,  on  nickel  is  usually  less  than 
several  hundred  Angstroms. 

The  capacitance  measured  for  Ni(OH)^  films  in  phosphate  solutions  appears 
to  mainly  involve  the  space  charge  layer.  Nevertheless,  the  double  layer 
would  be  expected  to  provide  some  contribution  to  the  capacitance.  Normally, 
the  space  charge  capacitance  is  smaller  than  the  double  layer  capacitance,  and 
so  dominates  the  overall  capacitance  (10).  However,  the  Mott-Schottky 
relationship  for  a  p-type  semi-conductor  predicts  that  the  space  charge 
capacitance  increases  as  the  potential  is  made  more  positive.  Accordingly, 


10 


work,  I  was  found  to  be  in  general  independent  of  potential  within 
experimental  error.  In  addition,  the  analysis  of  the  impedance  data  showed 
that  o  was  independent  of  potential  within  experimental  error  for  several  pH 
val ues . 

The  uncertainty  in  both  a  and  R  arises  from  several  sources.  The  only 
direct  contribution  is  from  the  calculation  of  these  parameters  which  are 
obtained  essentially  as  best-fit  values  and  accordingly,  have  a  given  error 
for  each  set  of  impedance  data. 

Various  other  factors  contribute  indirectly  to  the  uncertainty.  The  major 

one  is  the  difficulty  in  attaining  a  true  steady  state  for  the  growth  of  the 

films  and  this  may  have  been  at  least  partly  responsible  for  the  larger 

variation  in  o  with  potential  for  pH  values  of  7,  8,  and  10.  Although  the 

impedance  measurements  were  commenced  when  the  rate  of  change  of  the  dc 

current  was  negligible,  the  dc  "steady  state"  current  (typically  0.5  pA/cm2) 

generally  showed  some  decrease  over  extended  times.  Usually,  this  decrease 

o 

was  only  0. 1-0.2  pA/cm  but  it  did  indicate  an  inherent  problem  in  measuring 
the  impedance  of  a  film  at  various  potentials  when  the  film  thickness  is  not 
in  a  true  steady  state. 

A  further  factor  is  the  nature  and  composition  of  the  film  formed  at  each 
potential  and  pH.  It  was  found  in  previous  work  (4)  that  a  partially 
dehydrated  film  is  formed  on  nickel  in  a  Na?HP04  solution  at  pH=9.1.  This  film 
is  unlikely  to  be  uniform  in  composition,  but  probably  becomes  more  hydrated 
towards  the  fi lm/sol ution  interface.  Such  a  change  in  the  composition  would 
be  expected  to  affect  the  diffusion  characteristics  of  the  film,  thereby 
resulting  in  a  poorly-defined  diffusion  coefficient  For  the  mobile  species. 

In  addition,  the  degree  of  hydration  of  the  film  may  well  vary  quite 
significantly  over  the  pH  range  studied. 

The  final  aspect  to  be  considered  is  the  solubility  of  the  surface  film. 
Previously  published  thermodynamic  data  (11)  for  the  nickel/water  system  at 
25 °C  show  that  at  a  pH  slightly  greater  than  8,  Ni(0H)?  is  in  equilibrium  with 
dissolved  Ni2+  at  an  activity  of  lO-6  mol  kg'1.  Similarly,  at  a  pH  of 
approximately  12.5,  Ni(OH)^  and  Ni(OH)^  are  in  equilibrium.  However,  in  this 


Table  1 


Parameters  Derived  From  Impedance  Data 
For  0 -N i ( OH ) ^  Films  On  Single  Crystal  Nickel 


pH 

E/Y 

Cp/pF  cm  2 

R  /n  cm2 

P 

2 

o/n  cm  s 

7 

0.1 

16.8 

9573 

2859±143 

0.3 

19.5 

10956 

3486±174 

8 

0.1 

17.9 

6979 

3493H75 

0.3 

24.3 

3640 

5605±280 

9 

0.1 

16.4 

2946 

2827H41 

0.3 

24.6 

4593 

29121146 

10 

0.1 

18.1 

6948 

38591193 

0.3 

27.3 

2343 

57831289 

11 

0.1 

19.9 

5585 

37861189 

0.2 

24.7 

4831 

39151196 

0.3 

33.2 

4294 

39561198 

0.4 

40.6 

2833 

41161206 

12 

0.1 

20.7 

2715 

31541158 

0.3 

36.1 

2089 

30891154 

1/2 
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Values  of  the  various  circuit  elements  are  given  in  Table  1.  The  Warburg 
coefficient  is  independent  of  potential  within  experimental  error  for  pH 
values  of  9,  11,  and  12.  However  it  does  exhibit  a  steady  increase  with 
potential  at  each  pH  in  the  range  studied  (at  pH  values  of  7,  8,  and  10,  this 
increase  is  greater  than  the  experimental  error).  In  contrast,  except  for  pH 
values  of  7  and  10,  Rt  decreases  with  increasing  potential. 

The  capacitance  is  also  potential  dependent  which  is  consistent  with  data 
reported  by  Madou  and  McKubre  (5)  for  NifOH)^  films  in  KOH.  These  workers 
attributed  the  measured  capacitance  to  the  presence  of  a  space  charge  layer  in 
the  film,  although  the  values  appear  to  be  too  high  to  be  completely 
attributed  to  this  source.  Instead,  the  values  listed  in  Table  1  are  also 
characteristic  of  the  electrical  double  layer  suggesting  that  the  capacitance 
of  N i ( 0 H ) ^  films  in  phosphate  media  is  due  to  space  charge  in  the  solution 
(double  layer). 

The  values  derived  for  the  various  circuit  elements  were  used  to  calculate 
the  total  impedance  of  the  equivalent  circuit  used  in  the  analysis.  Typical 
Nyquist  plots  of  calculated  and  experimental  data  are  shown  in  Figures  1  to  4. 
The  two  sets  of  data  are  in  generally  good  agreement,  particularly  in  the  high 
and  low  frequency  regions.  At  intermediate  frequencies,  there  was,  in  some 
cases,  a  small  difference  between  the  two  sets  of  plotted  data.  This 
difference  might  indicate  the  need  for  an  additional  component  in  the 
equivalent  circuit,  but  more  likely  it  reflects  the  difficulty  of  representing 
a  region  of  joint  activation/diffusion  control  in  film  growth  by  "ideal" 
circuit  elements.  Nevertheless,  the  equivalent  circuit  used  appears  to 
adequately  represent  the  impedance  associated  with  film  formation  on  nickel  in 
phosphate  solutions. 

Di scussion 


Previous  studies  by  Macdonald  and  co-workers  (6,3)  have  shown  that,  for 
the  passive  film  on  nickel  in  a  Na2HP04  solution  (pH=9.1),  the  electric  field 
strength,  the  charge  transfer  coefficient  and  the  dc  steady  state  current, 

I^j.  are  independent  of  potential.  Under  these  cqnditions,  it  is  predicted 
from  the  point  defect  model  proposed  by  Chao  et  al_.  (6,3)  that  o  should  also 
be  independent  of  potential.  Over  the  pH  range  (7  to  12)  used  in  the  present 
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SECTION  2 


THE  ANALYSIS  OF  AC  IMPEDANCE  DATA  FOR  NICKEL  ELECTRODES 


The  growth  of  hydroxide  films  on  single  crystal  nickel  in  phosphate 
solutions*  (pH=7  to  12)  was  investigated  using  ac  impedance  and  ell ipsometric 
techniques, and  details  of  the  experimental  work  and  results  have  been 
presented  previously  (4).  A  subsequent  analysis  of  the  ac  impedance  data  has 
been  carried  out  and  the  results  are  as  follows. 

Analysis 


The  impedance  data  were  analyzed  in  terms  of  an  equivalent  circuit  model 
in  which  the  electrochemical  interface  can  be  represented  by  a  capacitance, 

C  ,  in  parallel  with  a  series  combination  of  a  charge  transfer  resistance,  Rt, 
and  a  Warburg  (diffusion)  impedance,  Z w.  The  latter  is  given  (1)  by 


Zw  *  (l-j)ouf(1/2)  (1) 

where  j=  -1,  uj=2irf,  f  is  the  frequency  of  the  ac  signal  and  a  is  the  Warburg 
coefficient.  The  values  of  the  circuit  elements  were  determined  by  jdecon- 
volution  of  the  experimental  data  in  a  series  of  steps.  Firstly,  the  contri¬ 
butions  of  the  uncompensated  solution  resistance  and  Cp  were  removed.  The 
solution  resistance  was  determined  from  the  high  frequency  intercept  of  the 
impedance  data  plotted  in  the  complex  plane,  while  Cp  was  obtained  from  a  Bode 
plot  (log  | Z |  versus  log  ,  where  Z  is  the  impedance).  The  real  component  of 
the  reduced  impedance  data  was  then  plotted  against  to”^2^,  and  a  linear 
regression  performed  to  yield  values  of  R^  and  from  the  intercept  and 
slope,  respectively. 


*  These  solutions  contained  varying  amounts  of  NaH2P04»H20,  Na2HP04,  and 
Na,P04*12H20  adjusted  to  give  a  total  phosphate  concentration  of  0.1  mol 
t-1  at  the  desired  pH. 
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SECTION  1 


INTRODUCTION 


Reactive  metals  such  as  aluminum,  chromium,  titanium,  and  zirconium,  all 
of  which  play  an  important  role  in  our  metals-based  civilization,  can  be  used 
in  practice  only  because  of  the  existence  of  a  thin  oxide  film  at  the  metal 
surface.  This  "passive"  film  serves  to  isolate  the  underlying  metal  from  the 
environment,  and  hence  is  responsible  for  the  kinetic  stability  of  metals  when 
exposed  to  oxidizing  atmospheres.  However,  passive  films  generally  are  not 
perfect  and  may  breakdown,  resulting  in  localized  corrosion,  (pittinq.  etc.). 
Breakdown  is  normally  induced  by  some  aggressive  species  (e„g„  Cl  )  in  the 
solution.  Even  if  breakdown  does  not  occur,  the  passive  film  will  continue  to 
grow  and  itself  may  slowly  dissolve  into  the  environment,  resulting  in  a 
steady  state  thickness  on  the  metal  surface. 

The  mechanism! s)  of  growth  of  passive  films  on  metal  surfaces  have  been 
under  extensive  investigation  since  the  middle  of  the  nineteenth  century. 

This  work  has  sought  to  provide  an  explanation  of  experimental  data  regarding 
the  kinetics  of  film  growth  and  the  conditions  under  which  breakdown  occurs. 

In  general ,  these  efforts  have  been  unsuccessful,  either  because  the  various 
models  proposed  are  unrealistically  simple  or  because  insufficient  data  have 
been  available  to  distinguish  between  them. 

In  this  final  report,  we  present  experimental  data  which  have  been 
obtained  and  describe  the  various  theoretical  developments  which  have  been 
achieved  throughout  the  program.  The  experimental  studies  include  impedance 
analysis  of  passive  films  formed  on  single  crystal  nickel  in  phosphate 
solutions,  as  well  as  our  recent  work  on  passivity  breakdown  of  single  crystal 
nickel  in  buffered  sodium  chloride  solutions.  The  theoretical  studies  involve 
the  further  development  of  the  point  defect  model  which  was  originally 
described  by  Chao,  Lin,  and  Macdonald  (1-3).  The  theoretical  work  is  far  from 
complete  but  we  believe  it  will  form  a  new  basis  for  the  design  of  future 
experiments  to  test  various  theortical  descriptions  of  passive  films  which  are 
currently  being  developed. 


1 


PREVIOUS  PAGE 
IS  BLANK 


LIST  OF  TABLES 


Table  No.  Page 

1  Parameters  derived  from  impedance  data  for  B-Ni(OH)-  films  4 

on  single  crystal  nickel. 

2  The  critical  pitting  potentials  ( Vc )  of  pure  nickel  (100)  22 

in  borate  buffer  solution  at  different  chloride 
concentrations  ([Cl-]). 

3  The  induction  time  (t)  of  pitting  corrosion_of  single  crystal  27 

nickel  (100)  in  borate  buffer  solution  ([Cl  ]  =  0.2  M)  at 
different  breakdown  overpotential  (&V). 

4  Experimental  and  theoretical  values  of  growth  of  passive  44 

film  solution. 

5  Experimental  and  theoretical  values  for  growth  of  nickel  48 

film  in  phosphate  solutions. 


Figure  No. 

Page 

36 

Influence  of  Kn  on  the  pitting  potential  as  a  f unci  ton 
of  %  molybdenum. 

92 

37 

Influence  of  the  concentration  of  cations  vacancies  on 
pitting  potential  and  incubation  times  as  a  function  of 
molybdenum  content.  aV=vc(xmq)  _vc^XMo=0^' 

93 

38 

Influence  of  complex  stoichiometry  (q)  on  the  pitting 
potential  and  incubation  time  as  function  of  molybdenum 
content.  AV=VC ( XMq ) - Vc ( XMq=0) . 

95 

t 


xi  i 


Page 

74 


27  Calculated  electric  field  strength  versus  film  thickness. 

Note  that  in  this  case  the  anion  and  ionized-anion  vacancy 
concentration  at  the  film/solution  interface  are  varied. 

In  all  cases,  the  deviation  from  neutrality  is  negligibly 
smal 1 . 


28  Calculated  field  strength  versus  thickness.  Note  that 

in  this  case,  the  concentration  of  positive  holes  Ap 
is  varied. 


29  Calculated  electric  field  strength  versus  film  thickness. 

In  this  case  the  concentration  of  electrons  An  is  varied. 


75 

76 


30  Calculated  electric  field  strength  versus  film  thickness.  77 

Note  that  in  this  case  the  anion  and  ionized-anion  vacancy 
concentration  at  the  film/solution  interface  are  varied 
and  the  concentration  of  cation  vacancy  at  interfaces  is 
null.  In  all  cases  the  deviation  from  electro-neutrality 
is  small . 


31  Calculated  electric  field  strength  versus  film  thickness.  78 

Note  that  in  this  case  the  anion  and  ionized-anion  vacancy 
concentration  at  the  metal  /film  interface  is  varied  and 

the  concentration  of  cation  vacancy  at  interfaces  is  null. 

In  all  cases  the  deviation  from  electro-neutrality  is  small. 

32  Calculated  electric  field  strength  versus  film  thickness.  79 

Note  that  in  this  case  the  cation  and  ioni zed-cation 

vacancies  concentrations  at  the  metal /film  interface  are 
varied;  and  the  concentration  of  ion  vacancy  is  null. 

In  all  cases  the  deviation  from  electro-neutrality  is 
negligibly  small. 


33  Calculated  electric  field  strrength  versus  film  thickness.  80 

Note  that  in  this  case  the  cation  and  ionized-cation 
vacancies  concentrations  at  the  film/solution  interfaces 
are  varied;  and  the  concentration  of  ion  vacancy  is  null. 

In  all  cases  the  deviation  from  the  electro-neutrality 
is  negligibly  small. 


34  Pitting  potentials  versus  %  molybdenum  in  18%-Cr  ferritic  90 

stainless  steel,  {•experimental  work  in  1  M  NaCl  solution 

at  258C  (24).  Theoretical  results  q=l;  ny  =  4xl020cm~3, 

Kn=10"4cm3  (O);  nv=4xl019cm"3,  Kn=10"17cm3  (D)0. 

35  Effect  of  Molybdenum  on  induction  time  for  pit-initiation  91 

according  to  the  solute  vacancy  interaction  model  (SVIM). 

q=l,  nv=4xl020cm“3,  «n=10~4cm3). 


xi 


SECTION  3 


FILM  BREAKDOWN  STUDIES 


General  Objective 


Pitting  corrosion  at  passivated  metals  is  a  direct  result  of  the  breakdown 
of  passive  films  on  the  surface.  For  several  decades,  the  subject  of  pitting 
corrosion  has  been  studied  extensively,  yet  none  of  the  theories  proposed  have 
been  able  to  account  for  all  of  the  experimental  results.  The  recently 
developed  point  defect  model  (1)  provides  a  quantitative  explanation  of 
passivity  breakdown  in  terms  of  electrochemical  phenomena.  However,  this 
model  lacks  experimental  support,  since  most  of  the  existing  data  for  pitting 
corrosion  were  not  obtained  from  carefully  controlled  experiments.  The 
purpose  of  this  study,  therefore,  was  to  test  the  point  defect  model  by 
collecting  passivity  breakdown  data  under  well-defined  experimental 
conditions. 

Previous  Work 


During  the  preliminary  portion  of  this  study  (13,14),  an  experimental 
apparatus  was  assembled,  a  combined  mechanical /electrolytic  polishing 
procedure  was  developed,  and  a  Teflon  shield  was  designed  in  order  to  reduce 
crevice  corrosion.  Several  film  breakdown  potentials  (Vb)  were  measured  in 
borate  buffer  solution  containing  chloride,  but  no  analysis  of  the  data  was 
carried  out  (4).  During  the  second  year  of  study,  small  changes  were  made  to 
the  experimental  system  and  included  using  a  PAR  173  potentiostat/gal vanostat 
with  a  PAR  179  digital  coulometer  module  instead  of  a  Wenking  66TS10 
potentiostat  and  a  PAR  379  coulometer.  The  present  most  sophisticated 
experimental  system  is  shown  schematically  in  Figure  5.  The  polishing  of  the 
single  crystal  nickel  electrode  follows  the  procedure  developed  in  the  first 
year,  except  that  a  final  0.05  pm  alumina  polishing  was  performed  in  order  to 
remove  any  film  formed  during  the  preceding  electro-polishing  step.  The  same 
Teflon  shield  was  utilized  to  reduce  crevice  corroson  but,  in  order  to  ensure 
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Figure  5.  Experimental  apparatus  for  passivity  breakdown 
studies  on  nickel  in  borate  buffer/NaCl  solutions. 
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a  constant  exposed  surface  area  (Figure  6),  the  use  of  microstop  was 
discontinued. 

A  major  change  made  in  the  experimental  procedure  was  the  growth  of  the 

passive  film  in  the  chloride-free  solution,  which  greatly  increased  the 

reproducibility  of  the  data,  and  resulted  in  better-defined  breakdown 

potentials.  The  sodium  chloride  solution  was  added  after  a  8-Ni(0H)9  film  had 

-3  c 

been  grown  to  a  constant  thickness  (charge  0=5.000x10  coulomb)  on  the  single 
crystal  nickel  (100)  in  the  chloride-free  borate  buffer  solution  at  100 
The  time  required  to  add  and  mix  the  sodium  chloride  solution  with  the  borate 
buffer  solution  was  short,  such  that  the  change  of  film  thickness  during  this 
period  was  negligible  (less  than  22).  Potential  sweep  (film  breakdown 
potential  measurements)  or  the  potential  step  (induction  time  measurements) 
type  of  measurements  were  performed  immediately  after  the  mixing  to  determine 
the  breakdown  potential  ( Vb )  or  induction  time  (t),  respectively.  Figures 
7(a)  and  8(a)  show  the  potential -time  profiles  for  these  two  experiments. 

Results  and  Discussion 


i)  Film  Breakdown  Potentials  (V.  ) 

b 

From  potential  sweep  experiments,  the  film  breakdown  potential  (Vb)  was 
measured  at  the  point  where  the  current  increases  sharply  on  the  anodic 
polarization  curve  as  shown  in  Figure  7(b).  In  the  earlier  work,  when  the 
passive  film  was  grown  in  the  chloride-containing  solution,  the  film  breakdown 
potential  was  difficult  both  to  determine  and  to  reproduce  satisfactorily  (4). 
In  the  more  recent  work,  however,  the  film  was  grown  in  the  chloride-free 
solution  and  the  breakdown  potential  was  better  defined  and  was  also 
reproducible  (Figure  9).  This  difference  in  the  behavior  of  Vb  indicates  that 
chloride  ions  modify  the  passive  film  during  the  growth  stage,  and 
consequently  influence  the  susceptibility  of  the  film  to  breakdown.  For 
example,  MacDougall  et  al .  (15)  claim  that  chloride  ions  penetrate  the  passive 
film  on  nickel,  but  this  conclusion  has  not  been  established  unequivocally 
because  of  the  difficulty  of  depth-profiling  thin  films  using  sputtering/ Auger 
Electron  Spectroscopy  and  Secondary  Ion  Mass  spectrometry  (SIMS). 
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Figure  6.  Shield  design  for  reducing  crevice  corrosion. 
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Figure  7.  (a)  Potential -time  profile  used  in  film  breakdown 

studies,  (b)  Schematic  current-potential  curve  illustrating 
film  breakdown. 
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Figure  8.  (a)  Potential -time  profile  used  to  determine 

induction  time  for  passivity  breakdown,  (b)  Schematic 
current-time  curve  illustrating  passivity  breakdown. 
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Figure  10.  Current-potential  curves  for  Ni(100)  in  buffered 
borate  solution  with  0.2N  chloride  at  different  sweep  rates. 
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Figure  11.  The  effect  of  sweep  rate  on  the  film  breakdown 
potential  of  Ni(100)  in  borate  buffer  solution  as  a  function 
of  chloride  concentration. 
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Table  2 

Critical  Pitting  Potentials  ( Vc )  for 
Nickel  (100)  in  Borate  Buffer  Solution  at 
Different  Chloride  Concentrations  ([Cl-]). 


[Cl-] 

(M) 

0.02 

0.05 

0.1 

0.2 

0.5 

vc  <VSCE)  0.620 

0.560 

0.412 

0.350 

0.350 

log  [cr]  (N) 


Figure  12a.  Relationship  between  the  critical 
potential  and  chloride  concentration  fro  Ni(lC 
borate  solution. 


■p 


Figure  12b.  Plot  of  the  logarithm  of  induction  time  (log  t) 
versus  the  inverse  of  overpotential  ( 1 / AV ) . 
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Ni  (100) ,  [cr]  =  0.2M 


B 


2.303RT 

aF 


where  a  is  the  fraction  of  the  applied  potential  contributing  to  the  potential 
drop  at  the  film/solution  interface,  and  F,  R,  and  T  have  their  usual  meaning. 
The  intercept  A  cannot  be  analyzed  because  no  data  are  available  for  the 
parameters,  u,  Jm  and  J*.  Furthermore,  the  available  data  for  B  do  not  permit 
an  accurate  estimate  of  a  at  this  time  although  the  data  given  in  Figure  12a 
indicate  that  the  value  is  of  the  order  of  0.27.  This  value  is  much  lower 
than  that  (0.878)  quoted  by  Chao  et  al^  (18)  for  polycrystall ine  nickel  in 
phosphate  solution  (pH  4.5).  It  is  important  to  note  that  both  the  solution, 
in  particular  the  pH,  and  the  crystalline  form  of  the  nickel  differ  in  the  two 
studies,  and  either  of  these  factors  could  well  account  for  the  above 
difference  to  a  large  extent. 

iii)  Induction  Time  (t). 

The  induction  time  (t)  of  pitting  corrosion  at  a  certain  breakdown 

overpotential  (aV  =  V  -  V  )  can  be  determined  from  the  current-time  curve 

“PP  c 

in  the  potential  step  experiment  (Figure  10(b)).  For  single  crystal  nickel 
(100)  surface  the  induction  time  was  measured  at  twelve  different 
overpotentials  in  the  borate  buffer  solution  containing  0.2  mole/liter  sodium 

_3 

chloride  while  other  parameters  such  as  film  thickness  (0  =  5.0x10  coulomb), 
film  formation  potential  (+100  mV),  temperature  (25°C),  etc.  were  kept 
constant.  The  induction  time  increased  very  rapidly  with  decreasing 
overpotential  (Table  3).  The  relationship  between  the  induction  time  and 
overpotential  can  be  presented  in  two  ways: 

a.  log  t  vs  1/ AV 

When  the  logarithm  of  the  induction  time  plotted  against  the  inverse 
of  the  overpotential  (Figure  12b),  a  linear  relationship  was  found  with  very 
good  correlation  coefficient  (0.98): 

log  t  =  0.61  +  0.071  ( 1/aV )  (4) 


The  intercept  at  1/aV=0  (or  AV->uo)  was  4  seconds. 


Table  3 


Induction  Times  (t)  Of  Nickel  Specimens  In 
([C1~]=0.2M)  At  Different  Overpotentials  (AV). 


Ni-(IOO) 


AV  (V) 

t  (sec) 

0.026 

2300 

0.031 

774 

0.035 

305 

0.041 

229 

0.051 

78 

0.06 

94 

0.071 

38 

0.091 

32 

0.  Ill 

22 

0.151 

17 

0.201 

9 

0.226 

4 
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The  model  developed  in  the  previous  sections  will  be  used  to  calculate 

ilm  properties  which  can  be  compared  with  available  experimental  data.  The 

txperimental  data  of  particular  interest  pertains  to  the  growth  of  passive 

:ilms  on  iron  in  borate  solutions  (1),  and  on  nickel  in  phosphate  solutions 

6).  To  obtain  film  properties  from  the  theoretical  first-order  point  defect 

Kidel ,  it  is  necessary  to  input  numerical  data  for  a  number  of  parameters. 

ialues  for  some  of  these  parameters  are  well-known.  However,  values  for 

jthers  where  chosen  by  adjusting  various  parameters  so  that  the  calculated 

:hickness  agreed  with  experiment  at  one  point;  normally  that  at  -0.4  V  after 

>ne  hour  of  polarization  for  which  the  film  thickness  L  =  L  .  Details  of 

o 

:hese  calculations  for  passive  films  formed  on  iron  and  nickel  are  presented 
in  the  following  sections. 

I.  Iron  (Fe-Feo0o) 

- - C~  J— 

The  parameters  that  are  required  for  applying  the  model  to  the  formation 
)f  oxide  films  on  iron  were  chosen  as  follows: 

i)  Experimental  Data  (1) 


a  =  0.728  (at  20*0 

from 

Eq.  (10) 

B  =  0.0144  V 

from 

Eq.  (10) 

pH  =  11.5 

from 

Eq.  (10) 

-?n  2 

D  =  2  x  10  cm  /s 

from 

Eq.  (5) 

Q  =  30  cm3/mol 

from 

Eq.  (4) 

Physical  Constants 

<NJ 

11 

O' 

from 

Eq.  (5) 

F  =  96487  coul /mol 

from 

Eq.  (6) 

R  =  8.314  J/mol °K 

from 

Eq.  (6) 

T  =  298  *K 

from 

Eq.  (6) 

41 


M 

n 


dX 

dt 


d2X 


V  = 


dt 


X 


n 


=  (t  +  nh)  and  increment  h  =  (t  -  t  The  Taylor  series  is  often  used 

10  n  n  - 1 

i  starting  a  solution  to  be  carried  on  by  some  other  method. 


A  more  accurate  method  of  extrapolation  uses  instead  of  a  straight  line,  a 
irve  represented  by  an  algebraic  equation  of  degree  higher  than  one.  This  is 
le  basis  for  what  is  known  as  the  Adams  method.  The  curve  used  in  this 
>thod  is  a  parabola  which  can  be  represented  by  an  algebraic  equation  of 
;cond  order  as  follows: 


Vl  '  xn  *  h 


♦  (1/21  («„  -  «„_!) 


'  a  similar  process  a  curve  of  fifth  degree 
luation  for  extrapolation: 


Xn+1  =  xn  +  h  Mn  *  <1/2^n-l  +<5/12>A  M 


can  be  used  with  the  resulting 
n_2  +(3/8)A3Mn_3  +  (251/720)a\_4 


ere 


AM  •  =  M 
n- 1  n 


n- 1 


A2Mn  „  =  AM  .  -  AM  « 
n-2  n-1  n-Z 

A3M  ,  =  a2M  0  -  a2M  , 
n-j  n-Z  n-J 


a4m 


n-4 


=  a3m 


n-3 


?se  are  the  successive  differences  in  the  quantity  M,  and  similar  succesive 
Fferences  for  x  are  also  defined.  When  the  thickness  of  the  film  is 
Iculated  from  the  model,  it  is  possible  to  then  compute  the  concentrations 
oxygen  ion  vacancies  at  each  point.  These  calculations  are  then  repeated 
•  a  different  number  of  discrete  units  to  simulate  growth  of  the  film.  The 
iber  of  slices  is  chosen  arbitrarily. 
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—  =  —  M  /exp  -(q2k2D/4)ti-exp  (qk/2)l.l  • 

ft  ",  2  1  -  j  L  J 


C(f/S)  + 


C(m/f )  exp  [(qk/2)L  ]  -  C(f/s) 


C(m/f)  exp  C ( qk/2 ) Lq ]  -  C(f/s) 


It  is  not  possible  to  obtain  an  analytical  solution  to  this  equation,  and 
hence  it  will  be  necessary  to  use  a  numerical  method  to  theoretically  evaluate 
the  kinetics  of  film  growth.  A  numerical  method  was  developed  for  an  Apple 
Il-e  microcomputer,  and  this  program  has  allowed  the  generation  of  information 
on  the  following:  i)  growth  of  the  film,  ii)  thickness  of  the  film  as  a 
function  of  applied  potential,  Vext,  and  iii)  concentration  of  ion  vacancies 
throughout  the  film.  In  generating  numerical  data  for  L  as  a  function  of  time 
and  applied  potential,  the  basic  problem  is  to  solve  the  first-order 
differential  equation 


subject  to  the  initial  condition  that,  at  t  =  0,  L  =  L  .  It  is  assumed  also 
that  the  function  f(L,t)  is  generally  continuous  and  single-valued,  and  has 
definite  derivatives.  The  numerical  analysis  is  carried  out  by  first 
assigning  discrete  values  to  the  properties  of  the  film  as  a  function  of 
distance  and  time,  followed  by  successive  extrapolation  to  determine  specific 
properties  at  each  point. 

Perhaps  the  most  obvious  way  to  extrapolate  from  one  step  to  another  is  to 
develop  a  solution  to  the  differential  equation  in  terms  of  a  Taylor  series. 
This  solution  has  the  form: 


r,  =  x  +  M  h  +  1/2  M  ,h2  +  1/6  M  „h3  + 
n+1  n  n  n  n 


where 


Y  =  epx 


C'jcoslqx)  +  C£sin(qx) 


where  p  =  -a/2,  q  =  >/b-a  -  (a2/4),  i  =  \/-l  and  q,  CJ  are  integral 
constants.  The  appropriate  combinations  of  these  equations  yield  the 
particular  solution  which  allows  the  computation  of  the  concentration  of 
vacancies  as  a  function  of  time  and  distance.  The  most  acceptable  solution, 
physically,  is  the  one  where  both  roots  m^  and  m^  are  equal.  Accordingly,  the 
general  solution  can  be  written  as  follows 

C(x,t)  =  (A  +  Bx)  exp  C(qk/2)x]  exp  [-(q2k2D/4)t] 

in  which  the  constants  A  and  8  are  obtained  from  the  initial  and  boundary 
conditions: 

C(0,0)  =  C(f/s) 


C(Lq,0)  =  C(m/f) 


where  Lq  is  the  initial  thickness  of  the  film.  Accordingly,  the  solution 
becomes 

C(x,t)  =  exp  [-(q2k2D/4)t]  Ic(f/s)  (] 


X  exp  r(qk/2)x] 


C(m/f )  exp  (qk/2)Lo  -  C{f/s) 


From  Equation  (4)  and  (5),  it  is  possible  to  obtain  information  about  the 
growth  of  the  film.  Substituting  Equations  (18)  and  (5)  into  (4)  we  can 
write : 
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2 

The  sign  of  f,  will  be  chosen  to  provide  a  realistic  physical  description  of 
the  growth  of  a  passive  film  on  a  metal  surface.  It  is  evident  that  for  any 
given  value  of  x,  as  the  film  grows,  the  concentration  of  VQ..  at  this  point 
will  decrease  with  time.  Accordingly,  the  minus  sign  corresponds  to  the 
physically  realistic  situation,  and  Eq.  (12)  becomes 

R(t)  =  A  exp  (-E2  t)  (14) 

where  £  is  a  real  number  and  A  is  an  integral  constant.  Three  particular 
solutions  of  Eq.  (13)  are  possible,  depending  upon  the  roots  of  the  quadratic 
equation  in  m: 

2 

m  +  am  +  b  =  0 

2 

where  a  =  -  qK  and  b  =  +£  /D.  Let  m^  and  m2  be  the  roots  of  the  equation 
above.  Accordingly, 

CASE  1. 

mj  ?  m 2*  and  and  m?  are  real  numbers 

Y  3  Cj  exp  (m^x)  +  C?  exp  (m2x)  (15) 

CASE  2. 

m^  =  nr^,  and  nij  and  m2  are  real  numbers 

Y  =  exp  (mjx)  [Cj  +  X'C^  (16) 


CASE  3. 

m^  ^  m2,  and  mj  and  m2  are  imaginary  numbers 
mj  =  p  +  iq 
m2  =  p  -  iq 
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nonsteady  state  are  presented,  in  which  the  electrical  field  strength  is 
considered  to  be  constant. 

Solution  of  Diffusion  Equations 

For  each  species  moving  with  the  passive  film,  Ficks'  second  law  can  be 
written  as 


3C 

at 


♦  D  2£.  -  Dqk  !£ 
3x 


i? 


where  D,  q,  K,  are  considered  to  be  independent  of  distance  and  time  (field 
strength  through  the  film  is  constant).  An  analytical  solution  is  developed 
by  using  the  separation  of  variables  technique.  This  yields  the  following 
equation: 

C(x,t)  =  Y(x)  R(t) 

and  substitution  with  Eq.  (6)  yields 


1  3R(t) 

R(t)  at 


D  32Y(x)  _  Dqk  3Y(x) 
Y(x)  3x2  Y ( x )  3x 


Because  the  two  sides  are  separately  pure  functions  of  R(t)  and  Y(x)  we  may 
equate  both  to  the  same  constant,  E2 


1  3R(t)  ^  r2 

TTTtT  =  *  c 


D 

Y  (x ) 


32Y(x) 


Dqk 
Y  ( x ) 


9Y(x) 

ax 


(13) 
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By  involving  the  assumption  of  equilibrium  at  the  two  interfaces,  it  is 
possible  to  specify  the  boundary  conditions  as  follows  (1): 

a)  Metal /film  interface. 


The  boundary  and  initial  conditions  for  the  equation  (6)  are 


C(m/f ) 


n,  2fKn  '  l2/<)4G2i 

—  exp  - - - 

n  RT 


(8) 


b)  Fi lm/sol ution  Interface 


n  AG.,, 

v  31 


C( f/s)  =  —  exp 
n 


-  4.606  pH 


where 

0f/s  =  aVext  +  8  pH  *  •f/s  (10) 

•m/f  =  n'a)Vext  "  6  pH  ‘  0f/s  +  0R  "  0f  (11) 

a  =  Constant 
B  =  Constant 
pH  =  pH  of  solution 

=  potential  drop  across  the  film 
Vext  =  Externa11y  applied  potential 
0R  =  Reference  potential 

0f/$  =  Initial  potential  at  film/solution  interface 

AG'j  =  Gibbs  energy  of  the  reaction:  Vg-  ♦  H.,0  2H+(aq)  +  Og 


In  the  following  section,  we  report  the  results  obtained  from  a  formal 

solution  (analytical  and  numerical)  of  the  diffusion  equation  for  the 
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(4) 


dL  a  . 

*  \  V 


where  n  is  the  molar  volume  per  cation,  n  is  Avogadro's  number,  and  is 

v  v0- 

the  flux  of  oxygen  vacancies,  Vg-  towards  the  metal/film  interface. 


The  flux  Jv  can  be  obtained  from  Pick’s  first  law  ,  as  follows: 
0" 


Jv0  =-D^+D^c 


(5) 


where  C  is  the  concentration  of  Vg--.  The  concentration,  C,  as  a  function  of 
distance  and  time  is  obtained  by  solving  Pick's  second  law 


3C 

at 


^  -  Dqk 
3xZ  Sx 


where  k  is  given  by 


(6) 


k  = 


_  P  d0 

ITT  (Tx 


(7) 


The  others  parameters  are  as  follows: 

D  =  Diffusivity  of  VQ- 
q  =  Number  of  charges  on  species,  Vg- 
R  =  Gas  constant 
T  =  Temperature 
P  =  Faraday  constant 

0  =  Electrical  potential  through  the  film 
x  =  Distance 
t  =  Time 
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Figure  13.  Schematic  of  the  working  and  reference 
electrode  in  a  typical  electrochemical  cell  (upper  figure) 
and  the  corresponding  potential-distance  relationships 
(lower  figure). 
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In  developing  the  first-order  model,  the  previously-employed  simplifying 
assumptions  were  involved  as  follows: 


i)  The  film  is  assumed  to  be  an  oxide  of  composition  MO 


X/2 


ii)  The  passive  film  contains  a  high  concentration  of  point  defects;  the 
species  present  are  V^'*  Vg...  e',  h ,  but  only  Vg..  is  responsible  for 
the  growth  of  the  film  (1). 


iii)  The  electric  field  strength  which  exists  within  the  film  is  assumed 
to  be  constant  and  buffered  by  mobile  charge  carriers. 


i v )  Interactions  between  point  defects  are  considered  to  be  negligible. 

v)  Equilibrium  exists  at  the  metal  film  and  film/solution  interfaces 
with  regard  to  the  concentrations  of  the  point  defects. 


However,  in  the  first  order  model  developed  here  it  is  not  assumed  that  the 
film  is  in  a  quasi-steady  state,  as  was  done  in  the  original  treatment  (1). 
Consequently,  the  movement  of  point  defects  is  described  exactly  by  Ficks' 
second  law  (Nernst-Pl anck  equation)  rather  than  approximately  by  the  Ficks 
first  law  expression. 


Figure  13  shows  schematically  the  relationship  between  the  working 
electrode  and  the  reference  electrode  in  an  electrochemical  cell,  and  the 
corresponding  potential /position  relationships.  According  to  assumption  (ii) 
above,  growth  proceeds  by  the  movement  of  oxygen  ion  vacancies  from  the 
metal/film  interface  to  the  film/solution  interface  or  of  oxide  ions  in  the 
reverse  direction.  Accordingly,  fi’m  growth  occurs  exclusively  at  the  metal/ 
film  interface,  with  this  interface  advancing  into  the  metal  phase.  This 
interface  is  considered  to  be  planar  and  perpendicular  to  the  direction  of 
film  growth.  The  mathematical  problem  presented  here  is  clearly  a  linear  one, 
since  growth  takes  place  in  one  direction  only.  Since  we  have  assumed  that 
movement  of  oxide  vacancies  only  is  responsible  for  growth  of  the  film,  the 
rate  of  film  thickening  (dl/dt)  may  be  written  as 


metal  surfaces.  The  theoretical  equations  were  simplified  involving  certain 
assumptions  (1),  as  a  first  step  in  the  development  of  the  model.  The 
theoretical  results  obtained  were  then  compared  with  experimental  data,  and 
the  agreement  was  sufficiently  good  to  warrant  further  development.  In 
formulating  this  model,  the  following  processes  are  considered  to  occur  at  the 
two  interfaces  during  growth  of  the  film. 


solution  film  metal 

2H  +  Oq  Vq"+  H2O  e'  +  Vg‘‘(^)  +  Mm  m 

MX+  +  V  .  ^  Mm  Xe*  +  Mm  ^  V  v.  +  M 

MX  V  ”  M  MX 

Any  general  model  is  characterized  by  the  following  features: 

i)  The  species  present  are  V  {metal  vacancies);  VQ-(oxygen  vacancies; 
e'  (electrons);  h  (holes). 

ii)  The  various  species  interact  with  one  another. 

iii)  Movement  of  these  species  within  the  passive  film  is  described  by  the 
Nernst-PI anck  equation. 

iv)  The  interfaces  exist  under  non-equilibrium  conditions. 

v)  The  fieM  strength  within  the  film  may  be  a  function  of  distance  and 
ti  me . 
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SECTION  4 


THEORETICAL  DEVELOPMENT  OF  POINT  DEFECT  MODEL 


In  parallel  with  the  experimental  studies  reported  in  the  previous 
sections,  the  development  has  continued  of  a  point  defect  model  for  the  growth 
and  breakdown  of  passive  films  on  metal  surfaces.  These  developments  stem  from 
the  initial  work  of  Chao  and  co-workers  (1-3),  and  involve  principally  the 
relaxation  of  certain  simplifying  assumptions  in  the  original  model  and  an 
extension  of  the  theories  to  take  into  account  the  semi-conductor  properties 
of  passive  films.  This  new  model  is  far  from  complete;  however,  the  basic 
concepts  have  been  sufficiently  well  developed  as  to  be  presented  in  detail  in 
this  final  report. 

The  theoretical  modeling  work  on  film  growth  has  been  carried  out  at  three 
levels  of  sophistication.  At  the  first  level  (Subtask  4.1)  only  one  mobile 
species  and  a  constant  electric  field  strength  is  considered,  whereas  in 
Subtask  4.2  six  charged  species  are  considered  with  the  film  in  a  quasi 
steady-state,  and  accordingly,  film  growth  is  not  included  directly  in  the 
mathematical  description.  This  model  was  developed  to  explore  the  dependence 
of  the  electrical  and  electronic  properties  of  passive  films  on  the 
concentrations  of  charge  carriers.  Finally  {Subtask  4.3)  the  general  point 
defect  model  is  extended  to  consider  the  effect  of  minor  alloying  elements  on 
passivity  breakdown.  This  latter  "solute-vacancy  interaction  model" 
represents  a  major  departure  from  traditional  theories  for  the  role(s)  played 
by  minor  alloying  elements  in  the  breakdown  of  passivity  at  metal  surfaces. 

4.1  FIRST  ORDER  POINT  DEFECT  MODEL 

Introduction 


Assumption 

A  model  based  on  the  movement  of  point  defects  in  an  electrostatic  field 
has  been  proposed  (1-3)  to  interpret  the  growth  behavior  of  passive  films  on 
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1}  Comparing  the  passivity  breakdown  behaviors  (V^)  of  the  films  grown 
in  chloride-free  and  chloride-containing  solutions  suggests  that 
chloride  ions  modify  the  properties  of  the  passive  film  such  that 
sharp  breakdown  is  not  observed. 

2)  The  film  breakdown  potential  is  found  to  increase  linearly  with  the 
sweep  rate  for  five  different  chloride  concentrations,  although  the 
slope  varies  with  the  chloride  concentration. 

3)  The  critical  pitting  potential  is  proportional  to  the  logarithm  of 
the  chloride  concentration  as 

Vc(VSCE}  =  0,24  “  0,22  loS[c1~] 

According  to  the  point  defect  model,  the  constant  a  for  single 
crystal  nickel  (100)  in  borate  buffer  solution  is  calculated  from 
equation  (2)  as  0.27. 

4)  Induction  time  is  obtained  for  chloride  concentration  of  0.2 
mole/liter  at  twelve  different  overpotentials.  The  relation  between 
the  induction  time  and  overpotential  qualitatively  agrees  with  the 
point  defect  model.  Quantitative  agreement,  e.g.  the  value  of 
constant  a,  cannot  be  determined  because  of  the  unknown  constant  S' 
in  the  equation  (5). 


b. 


log  t  vs  AV 


The  point  defect  model  derived  the  relationship  between  the  induction 
time  and  the  overpotential  as 


t  =  £' 


(5) 


where  T  was  the  contribution  of  transient  period,  x  was  the  number  of  charge 
on  nickel  cation,  and 


£'  = 


5/J*  u"x/2  ([Cl"])x/2  exp 


XFaAV 

~wr 


(6) 


where  J”,  u,  and  £  were  constants. 

Figure  12c  showed  the  plot  of  log  t  versus  AV.  At  large  AV  the  induction 
time  approached  the  limiting  value  of  t  which  was  estimated  at  4  seconds  in 
this  case.  Substituting  x=2,  t=4,  and  constants  R,T,  and  F,  equation  (5) 
became 


t  -  4  =  £'  [exp(39  x  AV  x  a)  -l]'1  (7) 

By  adjusting  the  values  of  a  and  £',  curves  like  the  one  on  Figure  12c  could  be 
drawn  which  had  fairly  good  agreement  with  the  experimental  data.  However, 
the  value  of  a  could  not  be  determined  unambiguously  because  the  constant  £ 
was  unknown. 

Summary 


During  the  second  year  of  this  research,  the  passivity  breakdown  of  single 
crystal  nickel  (100)  has  been  studied  under  carefully  controlled  experimental 
conditions.  The  results  may  be  summarized  as  follows: 
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Reasonable  Estimates 


iii) 


C(L  ,0)(m/f)=  1.42  x  1024  cm'3 
C(0,0)(f/s)=  1012  cm"3 
AG“  =  -  60,000  J/coul 


from  Eq.  (8) 
from  Eq.  (9) 
from  Eq.  (8) 


iv)  Values  chosen  such  that  agreement  exists  between  the  experimental  and 
theoretical  initial  thickness  as  a  function  of  applied  potential. 


Lo  =  20  A 
rf/s  =  0.  38  V 
0R  =  0.27  V 
0f  =  0.2  V 


from  Eq.  (11) 
from  Eq.  (11) 
from  Eq .  (11) 
from  Eq.  (11) 


The  parameters  that  were  adjusted  to  obtain  agreement  between  theory  and 
experiment  for  the  initial  thickness  of  the  film  given  by  Chao  et  al.  (1),  are 
as  follows:  Thickness  at  the  initial  potential  (-0.4  V)  after  one  hour  of 
polarization,  Lq,  initial  potential  at  the  film/solution  interface,  0^$; 
initial  concentration  of  ionized  vacancies  at  metal /film  interface, 
C(LQ,0)(m/f) ;  reference  potential,  0R;  and  the  potential  drop  through  the 
film,  0f.  Accordingly,  the  first  point  calculated  ‘'heoretically  (thickness 
versus  applied  potential)  corresponds  to  the  initial  point  in  the  plot 
obtained  from  the  experimental  data  quoted  by  Chao  et  al  (1).  The  subsequent 
values  of  thickness,  L,  for  a  polarization  time  of  one  hour  as  a  function  of 
applied  potential,  V  calculated  from  the  model  are  compared  with 
experimental  data  in  Table  4  and  in  Figure  14.  It  is  apparent  that  the 
theoretical  points  are  close  to  the  experimental  ones,  and  in  both  cases  the 
dependency  of  film  thickness  on  potential  is  linear. 


Figure  15  shows  the  growth  of  the  film  as  a  function  of  time.  The  film 
grows  at  the  metal /film  interface  upon  stepping  the  potential  from  -0.4  V  to 

o 

0.785  V  with  an  initial  thickness  of  L  =  20  A.  After  one  hour  of  growth  at  a 

O  o 

potential  of  0.785  V,  the  thickness  is  close  to  48  A.  The  film  starts  to  grow 
slowly,  then  more  rapidly,  and  finally  reaches  a  plateau  at  much  longer  times 
(not  shown).  This  behavior  is  in  agreement  with  experimental  observations 
(2). 
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The  concentration  of  oxygen  ion  vacancies  at  the  metal/film  interface 
during  the  growth  of  on  iron  under  the  same  conditions  shown  in  Figure 

15  is  shown  as  a  function  of  film  thickness  in  Figure  16.  The  increase  in  the 
interfacial  oxygen  ion  concentration  is  a  direct  consequence  of  Equation  (8) 
because  of  the  decrease  in  as  the  film  thickens.  Note  that  in  this 
first-order  model,  it  is  assumed  that  equilibrium  exists  at  both  interfaces 
with  respect  to  the  point  defects. 


Table  4 

Experimental  and  Theoretical  Values  for  the  Thickness  of  the  Passive  Film  On 
Iron  in  Borate  Buffer  Solution  After  1  hr  of  Polarization 
(See  Reference  (1)  for  a  discussion  of  the  Experimental  Data.) 


Experimental  Theoretical 


Potential 

(V) 

Thickness 

o 

(A) 

Potential 

(V) 

Thickness 

(A) 

* 

* 

★ 

* 

-0.4 

20 

-0.4050 

20.10 

-0.2 

24 

-0.2300 

24.11 

0.0 

28 

-0.0150 

29.11 

+0.2 

32 

+0.1960 

34.04 

+0.4 

36 

+0.4110 

39.07 

+0.6 

40 

+0.6515 

44.70 

+0.8 

44 

+0.7846 

47.84 

★ 


Values  adjusted  (-0.4050,  20.1)  to  experimental  first  point  (-0.4,  20) 
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Figure  15.  Theoretical  growth  of  the  passive  film  (Fe^) 

as  a  function  of  time,  for  iron  in  borate  buffer  solution. 
Condition  as  in  Figure  14. 


I0‘24x  Concentration  of  Oxygen  Ion  Vacancies 
at  the  Metal/Film  Interface  (cm*3) 


Figure  16.  Concentration  of  oxygen  ion  vacancies  at  the 
metal/film  interface  as  a  function  of  film  thickness 
according  to  the  first  order  point  defect  model.  The 
calculations  are  for  an  film  on  iron  in  borate  buffer 

solution  under  the  conditions  stated  in  Figure  14. 
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Experimental  and  theoretical  data  have  been  discussed  by  Macdonald  et  al_. 
(6)  for  the  growth  of  passive  films  on  nickel  over  a  wide  range  of  pH 
(1.5-12.6)  in  phosphate  solution.  The  results  show  that  the  solution  pH  has  a 
strong  influence  on  the  formation  and  growth  of  passive  films.  However,  it 
was  reported  that  in  basic  solutions  (pH:  9.1-12.6),  and  in  moderately  acidic 
solutions  (pH  =  4.5),  the  passive  films  possess  the  same  optical  constants, 
but  exhibit  different  growth  behaviors.  The  experimental  film  thickness 
values  were  plotted  against  time  in  accordance  with  the  logarithmic,  the 
inverse  logarithmic,  and  parabolic  laws,  for  different  values  of  applied 
potential (6) .  It  was  apparent  that  different  growth  laws  apply  during  the 
various  stages  of  growth. 

Using  the  rate  law  developed  from  the  point  defect  model,  and  by  comparing 
the  experimental  data  presented  in  reference  (6)  with  this  rate  law,  it  is 
possible  to  fit  the  parameters  used  in  the  theoretical  model  in  order  to 
obtain  an  estimate  of  the  film  thickness  at  the  initial  time  (1  minute  of 
polarization).  Table  5  summarizes  values  for  film  thickness  at  different 
times,  as  computed  from  the  experimental  data  reported  by  Macdonald  et  al  (6), 
as  well  as  theoretical  values  obtained  from  the  first  order  point  defect 
model.  The  values  presented  in  Table  5  correspond  to  passivation  of  nickel  in 
a  0.05M  NaH2P04  (pH  =  4.5)  solution  at  25 °C  and  at  an  applied  voltage  of 
1.24  V  (vs  SCE). 

In  order  to  adjust  the  initial  thickness  (i.e.,  that  at  1  min.)  calculated 
from  the  model,  so  that  it  coincides  with  the  experimental  data,  the  values 
chosen  for  the  parameters  are  as  follows: 

Electric  field  strength  (e)  =  1.4  x  108  V/cm 
Initial  thickness  ( L0 )  =  16  x  10"8  cm 
Coefficient  of  diffusion  (D)  =  1.3  x  10~17  cm2/s 
Concentration  of  oxygen  ion  vacancies: 


w 


(i)  At  film/ solution  interface,  C( f / s )  =  1012  cm”3 

(ii)  At  metal/film  interface,  C(m/f)  =  0.74  x  1022  cm  3 


Table  5 

Experimental  and  Theortical  Values  for  Growth  of 
Nickel  Oxide  Film  on  Nickel  Phosphate  Solution 
{V  =  1.24  V  vs  SCE,  pH  =  4.85,  T  =  25*0 


From  Reference  (6) 

a 

Thickness  (A)  t(min) 

From  first  order 
• 

Thickness  (A) 

Point  Defect  Model 

t(min) 

16.3 

1 

16.3 

1 

18.5 

3 

16.9 

3 

25.0 

10 

19.2 

10 

27. 5 

15 

20.8 

15 

29.5 

20 

22.6 

20 

33.0 

30 

29.1 

30 

All  of  these  values  are  considered  to  be  physically  reasonable.  Also,  the 
electrical  field  strength  is  of  the  order  of  that  required  for  electron/hole 
pair  generation  and  subsequent  field  buffering,  as  assumed  in  the  point  defect 
model,  and  the  initial  thickness  reasonably  approximates  that  for  an 
air-formed  film.  The  concentration  of  oxygen  ion  vacancies  at  the 
film/solution  interface  is  very  low,  as  expected,  because  of  their  reaction 
with  water  molecules.  The  concentration  of  oxygen  ion  vacancies  at  the 
metal /film  interface  corresponds  to  about  25%  of  the  oxygen  ion  sites  being 
vacant;  a  requirement  that  also  seems  to  be  eminently  reasonable.  Figure  17 
shows  experimental  and  theoretical  plots  for  the  growth  of  the  passive  film  as 
a  function  of  time.  Although  significant  differences  exist  between  both  data 
sets,  the  form  of  the  change  in  film  thickness  with  time  is  the  same  in  both 
cases.  It  is  assumed  that  better  agreement  could  have  been  obtained  by 


\  ' 

wr’. 
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Time,  min. 

Figure  17.  Theoretical  and  experimental  data  for  the  growth 
of  the  passive  film  on  nickel  in  0.1  N  Na2H2P04  (pH  =  4.85) 

solution  at  25'C  and  at  an  applied  voltage  of  1.24  (V  vs 
SCE). 


further  refinement  of  some  of  the  parameters  in  the  model,  but  this  would  have 
proven  little  of  physical  significance. 

Because  some  of  the  parameters  used  in  calculating  the  film  thickness  from 
the  first  order  point  defect  model  are  not  well-known,  it  is  necessary  to 
carry  out  a  sensitivity  study  of  the  model  in  order  to  determine  how  the 
calculated  film  thickness  responds  to  each  of  the  parameters  of  interest. 

This  is  done  by  varying  each  of  the  parameters  in  turn  over  a  wide  range. 

Sensitivity  Study 

The  sensitivity  study  was  carried  out  using  the  mathematical  model 
described  in  the  previous  section.  The  values  of  some  of  the  parameters  used 
for  calculating  the  properties  of  the  passive  film  formed  on  iron  in  borate 
solution  were  varied  systematically,  and  the  rate  of  growth  of  the  film  and 
changes  in  the  concentration  of  oxygen  ion  vacancies  were  observed.  The 
parameters  chosen  for  variation  in  this  analysis  were  the  oxygen  ion  vacancy 
concentration  at  the  metal  film  interface,  and  the  vacancy  diffusivity.  These 
parameters  were  varied  over  the  ranges  of  1  x  1022  cm-3  to  5  x  1023  cm-3  and 
1.3  x  10  24  cm2/s  to  1.3  x  10  16  cm2/s,  respectively.  The  values  for  the 
other  parameters  remained  the  same  as  those  in  the  previous  section  for 
Ni(OH)^  films  on  nickel. 

Figure  18  shows  the  sensitivity  of  the  rate  of  growth  of  the  film  to 
variations  in  the  diffusion  coefficient.  Clearly,  the  diffusion  coefficient 
can  have  a  strong  effect  on  the  growth  kinetics  of  the  passive  film  on  nickel, 
at  least  in  the  case  considered.  Also,  as  shown  in  Figure  19,  a  range  exists 
in  the  value  of  the  diffusion  coefficient  over  which  this  parameter  has  a 
profound  effect  on  the  kinetics  of  film  growth.  This  range  appears  to  include 
the  likely  values  for  the  diffusivity  of  oxygen  ion  vacancies,  although  some 
data  indicate  that  D  is  of  the  order  of  10"21  cm2/s  (1). 

A  better  understanding  of  the  behavior  observed  in  Figure  19  is  obtained 
if  the  structure  of  Equation  (19)  is  studied,  which  relates  the  growth  of  the 
film  to  the  diffusion  coefficient.  Equation  (19)  can  be  written  as  follows 
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Figure  18.  Variation  of  film  thickness  and  oxygen  ion 
vacancy  concentration  as  a  function  of  time,  for  different 
values  of  the  vacancy  diffusion  coefficient. 


Diffusion  Coefficient ,cme/s 


Figure  19.  Calculated  thickness  of  the  passive  film  on 
nickel  as  a  function  of  the  diffusion  coefficient. 
Polarization  time  =  30  min,  other  parameters  as  indicated  in 
the  text. 


^  =  AD (8  +  CL)  exp(-aDt)  exp(0L) 


^here,  for  the  present  purposes.  A,  B,  C,  a,  and  6  will  be  considered  to  be 
:onstants.  Equation  (19)  predicts  the  following  behaviors: 

i)  If  D  — >  0 — >  dL/dt  — >0  (i.e.,  the  film  doesn't  grow,  if  L 

remains  finite) 

ii)  If  D  — ^  r<  dL/dt  — >0  (again,  the  film  doesn't  grow,  if  L 

remains  finite) 

iii)  There  exists  a  value  for  D  such  that  dL/dt  goes  to  a  maximum.  This 
value  is  identified  in  Figure  19. 

>le  have  yet  to  fully  appreciate  the  physical  implications  (if  any)  of  (ii)  and 
(iii)  above,  but  that  indicated  by  (i)  is  quite  obvious.  It  may  well  be  that 
(ii)  is  a  mathematical  artifact;  a  possibility  that  should  be  considered  in 
future  work.  Figure  20  shows  the  thickness  as  a  function  of  time  for 
Jifferent  values  of  initial  concentrations  of  oxygen  vacancies  at  the 
netal/film  interfaces.  Although  this  parameter  is  strictly  not  an  independent 
/ariable  since  it  depends  upon  the  applied  voltage  (see  Equations  (8)  and 
(ID)  it  is  apparent  from  the  data  in  Figure  20  that  the  rate  of  film  growth 
is  extraordinarily  sensitive  to  variations  in  the  concentration  of  oxygen  ion 
/acancies  at  the  metal/film  interface. 

<inetics  of  Film  Growth 


Over  the  past  several  decades,  a  great  deal  of  effort  has  been  invested  to 
jxplain  the  kinetics  in  the  growth  of  anodic  passive  films.  A  number  of 
jmpirical-theoretical  rate  laws  have  been  proposed;  the  best  known  being  the 
'logarithim  law"  and  the  "inverse  logarithim  law".  In  this  section,  a 
:omparison  is  presented  of  the  two  rate  laws  in  terms  of  the  mathematical  and 
)hysical  structure  of  the  first  order  point  defect  model  presented  earlier  in 
this  report. 
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general  analysis  of  Equation  (18)  allows  the  prediction  of  the  form  of 
ate  of  growth  of  the  film  as  a  function  of  time.  From  previous  analysis, 
ion  (18)  can  be  written  as  follows: 


C(x,t)  =  (A  +  Bx)  exp(-Bt)  exp(ax)  (21) 

s  A,  B,  8,  and  a  are  considered  to  be  constants  for  the  present  purpose, 
rdingly,  the  growth  of  the  film  is  expressed  as 

dL/dt  =  (A  +  BL)  exp(-o't)  exp(BL)  (22) 


orm  of  the  equation  shows  that 

i)  for  t->+O0,  dL/dt  — ^  +  0 

i)  t— >+  0,  dL/dt — ^constant 

P  P 

i)  d  L/dt  =  0,  inflexion  point 

eneral  shape  of  the  curve  L  =  f(t)  can  be  presented  by  the  plot  shown  in 
e  21.  The  form  of  the  rate  equation  shown  in  this  figure  will  now  be 
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he  algorithim  used  in  these  calculations  was  developed  for  and  executed 
Apple  Ile^  microcomputer.  The  algorithim  was  executed  for  both  sets  of 
ary  conditions  discussed  above,  and  the  resultant  electric  field  strength 
altage,  A 0  (where  A0  =  [eJ(x)  -  e°]Ax)  are  plotted  in  Figures  24  and  25, 
stively.  This  latter  quantity  is  a  measure  of  the  deviation  of  the 
tial  from  the  constant  field  approximation,  and  hence  is  a  direct 
ator  of  the  importance  of  space  charge.  From  Figures  24  and  25  it  is 
ant  that  space  charge  is  far  less  important  when  the  electrical 
ality  constraint  is  invoked  than  when  it  isn't.  This  same  behavior  was 
by  many  sets  of  data,  the  results  of  which  are  not  displayed  here  (see 
section) . 


Jemark  of  Apple  Computer,  Inc.,  Cupertino,  California 
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-2  CD(0)  -  CD+(0)  +  2  CA(L)  +  Ca_(L)  -An  +  Ap  '  0 

ere  L  represents  the  metal/film  interface  and  0  the  film/solution  interface. 

As  far  as  the  boundary  conditions  are  concerned,  the  concentrations  of 
ch  of  the  species  at  the  two  interfaces  must  be  specified.  Accordingly,  the 
undary  conditions  are  simply  stated  as  follows: 

24  -3 

An:  density  of  electrons  across  the  film  3  10  m 

24  -3 

Ap:  density  of  holes  across  the  film  =10  m 

CD+(0):  density  of  Cp+  at  f/s  interface  =  1024  m  3 

24  -3 

CD (0 ) :  density  of  CD  at  f/s  interface  =10  m 

CA_(0):  density  of  CA~  at  f/s  interface  =  1014  m  3 

CA (0 ) :  density  of  CA  at  f/s  interface  =  lCr  m  J 

14  -3 

Cd+(L):  density  of  CQ+  at  m/f  interface  =10  m 

Cd(L):  density  of  CD  at  m/f  interface  =  1014  m  3 

24  -3 

C.  (L):  density  of  C.-  at  m/f  interface  =10  m 

A~  A  24  -3 

Ca(L):  density  of  CA  at  m/f  interface  =10  m 

;se  values  have  been  chosen  such  that  electrical  neutrality  is  obeyed. 

A  second  set  of  boundary  conditions  was  chosen  without  regard  for 
;ctrical  neutrality.  This  was  done  to  explore  the  commonly-involved 
:uation  in  calculations  of  this  type  of  not  requiring  the  film  as  a  whole  to 
electrically  neutral.  The  boundary  conditions  chosen  in  this  case  are  as 
lows : 

An  =  1024  rrf 3 
Ap  =  1025  m" 3 
D+(0)  =  1024  m"3 
CD(0)  =  1024  nf3 
A_(0)  =  1024  m"3 
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y  number  of  charges  for  specie  CQ+  =  1 

2‘.  number  of  charges  for  specie  =  2 

number  of  charges  for  specie  C^-  =  -1 

4:  number  of  charges  for  specie  CA  =  -2 

:  initial  thickness  =  1A 
o  . 

p.-  Final  thickness  =  5A 

):  Potential  at  film/solution  interface  =  0.01  V 
):  Potential  at  metal /film  interface  =  0.05  V 

that  the  thickness  of  the  film  has  been  set  at  5A,  which  represents  the 
nsion  expected  for  the  very  thin  films  which  form  on  metals  and  alloys 
i  as  nickel*,  iron  and  stainless  steels  under  some  conditions. 

In  addition  to  the  normal  boundary  conditions,  the  film  must  be 
:trically  neutral,  i.e.. 


L  L  L 


luse  only  the  steady  state  condition  is  considered,  the  concentrations  are 
'  functions  of  x,  so  that  the  partial  differentials  become  normal 
erentials  and  it  is  possible  to  write  Equation  (35)  as: 

VCD(U  '  CD(0)]  +  qD+CCD+(L)  '  CD+(0)]  +  qACCA(L)  “  CA(0):i 
+  qA_[CA_(0)  -  CA_(0)]  -  eno  +  epQ  %  0 

tese  values  were  deduced  from  the  experimental  data  obtained  for 

issive  films  on  Ni  (18),  A0/ax  =  10^  V/cm. 
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.  -  JJ 


-^'tep0  exp(-B0)  -  enQ  exp(B0)  *  exp  t-Q^-j^-)  *3 


+  A-  exp[-q?k( — )x] 
9x 


+  A3  exp[-q3k{ — )x]  +  A4exp[-q4k( — )x] 
3  x  9  x 


qlBl  +  q2B2  +  q3B3  +  Q4B4 


Equation  (28)  can  be  expressed  in  general  form  as 


3 ct  /30 


— *■  =  fl — ,0,x 

3x^  \3x 


It  is  possible  to  construct  an  algorithim  to  solve  for  0,  in  which  the  initial 
filed  strength  (e°)  is  considered  to  be  constant.  The  value  obtained  is  then 
used  in  Equation  (33)  to  determine  a  new  value  for  e  (i.e.  e(x)),  and  so  forth 
for  j  iterations.  The  and  e^”1^  values  are  compared  and  if  the  difference 
is  less  than  a  preselected  tolerance,  the  value  e J (x )  is  taken  as  representing 
the  electric  field  strength  within  the  film  at  point  x. 

In  order  to  employ  the  algorithim  outlined  above,  we  must  choose  values 
for  the  various  constants  used  in  the  equations  and  specified  in  the  boundary 
conditions.  Also  the  tolerance  that  will  be  used  in  the  numerical  model  must 
be  specified. 

The  constants  of  interest  are  as  follows: 


F:  Faraday  constant  =  96487  coul/mol 

R:  general  constant  of  gases  =  8.314  Joule/mol #k 

T:  temperature  =  298 #K 

•*> :  dielectric  constant  for  oxide  iron  film  =  3.9 
F, :  dielectric  constant  of  vacuum  =  8.85419  X  10  ] 

_1Q 

e:  electron  charge  -  1.6  X  10  coulombs 

-23 

K:  Boltzman  constant  =  1.38  x  10  Joule/°kmol 
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(27) 


V20  =  (1/C -6)  Cq0+*CD+(x)  +  qD-C0(x)  +  qA‘CA(x)  +  qft-  -CA~(x) 
-  e.n  +  e.p] 


here  C  is  the  dielectric  constant  of  the  film.  At  steady-state,  Fick's 
econd  law  for  the  diffusion  and  electromigration  of  charged  species  can  be 
ritten  as  follows 


32C.  3C. 

0  =  D i  — J  -  D.^k  — - 
1  3x^  1  3x 


(28-31) 


here  i  =  1, 2,3,4  (i.e.  for  species  Cp,  Cp+,  Cft,  C A~).  Accordingly,  the 
lathematical  problem  of  solving  for  the  seven  unkowns  ( ,  C 2,  C3,  C4,  n, 

,  e)  requires  seven  independent  equations  given  by  Equations  (24-31).  It  is 
mportant  to  note  that  Equations  (27-31)  involve  coupled  non-linear  second 
rder  differential  equations.  These  must  be  solved  numerically. 

As  a  first  step,  consider  that  e  =  -d0/dx  =  constant.  Therefore,  from 
quations  (28-31)  it  is  possible  to  obtain  a  solution  for  the  concentration  of 
ach  species  of  the  form 

Cj  =  A,,  expfq^k  ex)  +  B.  (32) 

here  A.  and  are  integral  constants  and  are  functions  of  the  boundary 
onditions  for  each  species  considered. 

Substitution  of  Equation  (32)  into  Equation  (27)  permits  the  derivation  of 
non-linear  second  order  differential  equation  which  may  be  solved 
umerically  for  the  potential  0 
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The  species  considered  are 


concentration 

species 

charge 

CD 

donors  (V  Y < ) 

MX 

-  e 

V 

ionized  donors  (v^(X-i)^ 

-(x-De 

CA 

acceptors  (VQ--) 

+2e 

CA- 

ionized  acceptors  ( VQ") 

+e 

n 

electron  (e) 

-e 

P 

holes  (h) 

+e 

Numerical  Analysis 


Consider  that  6  species  are  involved  in  various  reactions  that  take  place 
at  the  film/solution  and  metal/film  interfaces.  The  species  correspond  to 
(^(acceptors),  CA-(ionized  acceptors),  CD(donors),  CD+(ionized  donors), 
electron,  and  holes.  It  is  supposed  that  the  concentrations  of  the  ionized 
and  unionized  vacancies  (CA,  CA~,  Cg,  CQ+)  follow  Fick's  second  law  of 
diffusion,  and  that  electron  and  holes  exist  under  equilibrium  conditions. 
Hence,  their  concentrations  are  given  to  first  approximation  by  the  Boltzman 
expressions  (28) 


n  =  nQ  exp  (+60) 

(24) 

P  =  PQ  exp  (-B0) 

(25) 

where  n  is  the  concentration  of  electrons  through  the  film,  P  is  the 
concentration  of  holes  through  the  film;  nQ,  P  ,  and  8  are  constants;  and  0  is 
the  potential  at  any  given  point  in  the  film.  The  electrical  potential  in  the 
oxide  is  governed  by  the  density  of  free  charge  (i.e.  the  "space  charge") 
according  to  Poisson's  equation: 

v20  =  -  4  <26) 


or 
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The  species  considered  are  donors  (V  x<),  ionized  donors  (V  (x-1)'^’  accePtors 
v  M  .M 

( V q— ) >  ionized  acceptors  (Vg..),  electrons  (e1)  and  holes  (h).  These 

calculations  permit  the  estimation  of  the  field  strength  through  the  film,  and 

to  determine  how  the  field  strength  is  affected  by  changing  the  concentrations 

of  various  species,  i.e.  by  varying  the  boundary  conditions.  Also,  the 

results  allow  the  determination  of  the  conditions  under  which  space  charge 

within  the  film  is  important,  and  hence  to  specify  the  conditions  under  which 

the  constant  field  strength  hypothesis  might  be  valid.  It  should  be  mentioned 

that  the  results  obtained  in  this  sub-task  will  be  applicable  only  to 

steady-state  films,  and  cannot  be  used  to  examine  the  effect  of  space  charge 

on  growing  oxide  films. 


Since  electrons  and  holes  are  obviously  present  in  anodic  films,  and  are 
intimately  associated  with  the  presence  of  donor  and  acceptor  impurities,  it 
is  fitting  to  enquire  into  the  roles  that  all  of  these  species  may  play  in 
determining  the  electronic  and  defect  structure.  The  model  adopted  for  this 
purpose  can  be  represented  as  follows: 

solution  film  metal 

chemical  reaction  Mm  ^  v  +  ^  m  +  Mm  + 

chemical  reaction  Vg»+  ^O^Og  +  2H  m  +  (x/2)  Vg  +  Xe 

electron  reaction  Mm  V  x_j  +  M*+  +  e1  m  +  V '  (x-D^^m  + 

electron  reaction  Vg-  +  H^O  ^  Og  +  2H  +  e'  m^M^  +  (x/2)  Vg  +  x/2  e 

x  =  L 
t  *  t 


x  =  L 
*  o 

t  *  0 
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in  the  model,  including  anion  vacancies  (VQ..),  cation  vacancies  (V  ^,), 

electron  { e ' )  and  holes  (h).  Additionally,  it  is  proposed  that  the  oxide 
possesses  the  normal  electronic  band  structure  which  characterizes  classical 
semiconducting  materials,  but  that  the  anion  vacancies  and  cation  vacancies 
act  as  electron  acceptors  and  donors,  respectively.  Thus,  in  terms  of  energy, 
these  species  are  assumed  to  occupy  positions  within  the  band  gap,  and  because 
the  concentrations  of  the  vacancies  are  very  high  (10  cm  vs.  10  cm  for 
dopants  in  a  classical  semiconductor)  the  electonic  properties  of  the  film  are 
expected  to  be  dominated  by  the  defect  structure.  The  picture  is  further 
complicated  by  the  fact  that  the  vacancies  are  not  uniformily  distributed 
across  the  film.  Thus,  in  our  view,  a  passive  film  is  best  described  as  a 
highly-doped  n-type/p-type  semiconductor  with  the  material  close  to  the 
metal/film  interface  being,  on  balance,  p-type  in  character  whereas  that  at 
the  film/ solution  interface  is  thought  to  be  n-type  in  character.  Clearly, 
the  model  that  emerges  from  these  considerations  is  a  complex  one,  but 
nevertheless  is  one  which  can  be  developed  mathematically  using  techniques 
which  have  been  employed  in  semi-conductor  physics. 

4.2  ROLE  OF  SPACE  CHARGE  IN  PASSIVE  FILMS 

Introduction 


Space  charge  in  metal  oxide  films  has  long  been  recognized  by  some  authors 
(24,25),  but  not  others  (26,27).  This  controversy  has  yet  to  be  resolved, 
principally  because  of  the  lack  of  a  direct  experimental  method  for 
determining  the  distribution  of  charge  in  very  thin  films. 

In  this  section,  the  distribution  of  charge  within  an  anodic  oxide  film 
using  the  point  defect  model  under  steady-state  conditions  is  calculated. 

These  calculations  differ  from  those  previously  reported,  in  that 
a  much  larger  number  of  charged  species  are  considered.  A  model  involving  six 
species  under  steady-state  conditions  in  which  the  electric  field  strength  is 
considered  to  be  a  function  of  distance  is  presented.  The  distribution  of 
each  species,  and  hence  the  space  charge,  is  obtained  by  solving  Fick's  second 
law  equation  for  each  species  coupled  to  the  Poisson  -  Boltzmann  equation. 

62 


(iii)  The  model  also  explains  why,  In  some  cases,  the  semi  empirical 
logarithmic  and/or  inverse  logarithmic  laws  can  be  fitted 
satisfactorily  to  the  same  experimental  data. 

From  the  sensitivity  study,  the  model  predicts  that  the  kinetics  of  growth 
of  the  anodic  film  on  nickel  is  a  strong  function  of  the  diffusion  coefficient 
of  the  oxygen  ion  vacancy.  It  is  well  known  that,  in  ionic  crystals  with 
impurities,  the  diffusion  coefficient  of  defects  (including  vacancies)  is  a 
function  of  both  temperature  and  concentration  (20,21).  Since  the 
concentration  of  vacancies  changes  throughout  the  film  and  is  a  sensitive 
function  of  the  intitfal  boundary  condition,  it  is  likely  that  kinetics  of  the 
film  growth  will  also  depend  upon  those  same  factors  which  determine  the 
vacancy  concentration  profiles.  This  may  explain,  at  least  in  part,  why  poor 
reproducibility  is  often  observed  in  anodic  film  growth  studies.  That  is, 
experimentally  it  is  extremely  difficult  to  generate  films  that  have 
reproducible  concentration  profiles  and  hence  vacancy  diffusion  characteristics. 
Very  few  experimental  data  exist  (22,23)  which  allow  the  establishment  of  a 
relationship  between  the  diffusion  coefficient  and  vacancy  concentration,  and 
hence  between  diffusivity  and  growth  of  the  film.  It  is  believed  that  these 
relationships  will  have  to  be  developed  if  more  accurate  descriptions  of 
passive  films  are  to  be  devised.  One  experimental  technique  which  shows 
considerable  promise  for  studying  the  motion  of  defects  in  passive  films  is 
the  AC  impedance  method  described  earlier  in  this  report  and  discussed  at 
length  by  Chao,  Lin  and  Macdonald  (1). 

Finally,  it  is  important  to  recognize  the  limitations  of  the  first  order 
point  defect  model  developed  in  this  section.  Because  only  a  single  mobile 
species  is  considered,  it  is  not  possible  to  account  for  space  charge  effects 
in  a  realistic  manner.  Also,  the  model  does  not  consider  the  electronic 
structure  of  the  film,  so  that  the  semi-conducting  properties  of  many  passive 
films  cannot  be  accounted  for.  In  order  to  include  phenomena  such  as  variable 
field  strength,  space  charge,  dissolution  of  the  film,  and  optical, 
electrical,  and  thermal  excitation  of  electrons,  it  is  necessary  to  develop  a 
more  advanced  model  which  combines  the  physical  and  electronic  properties  of 
the  film.  In  doing  so,  all  charged  species  must  be  included 
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Thickness 


( i )  t  ->  +  0 


1  nt  ■ 1  ^  -  oo 


L — y+  0 


(ii)  t — >.<  exp  (C/D)  lnt — >  <  oo  L“4  +  oo 

(ill)  t  >  exp  ( C/D )  1  nt  — ^  >  oo  L  — ^  -  oo 

(iv)  t  —yl  lnt — 0  L — >  1/C 

(v)  t — }*4  oo  lnt — +  co  L - ^  -  0 

These  conditions  are  described  by  the  general  function  plotted  in  Figure  23. 
Note  that  this  function  is  discontinuous  at  t=exp(C/D).  It  is  suggested  that 
the  only  domain  of  physical  significance  is  that  for  t<exp(C/D);  in  which  case 
the  growth  law  strongly  resembles  that  predicted  for  short  times  ( t<t ' )  by 
Equation  (22)  (see  Figure  21). 

Discussion 

The  results  obtained  from  the  first  order  point  defect  model  developed  in 
this  work  are  sufficiently  encouraging  to  warrant  that  the  model  be  further 
explored  (see  later).  Specifically,  the  simple  first  order  variant  has 
displaced  the  following  features: 

(i)  The  calculated  film  thickness  as  a  function  of  Vgxt  follows  a 
relationship  which  is  in  good  agreement  with  that  obtained  from 
experimental  data  for  iron  in  borate  solution. 

(ii)  The  rate  of  growth  of  the  passive  film  on  nickel  in  phosphate 

solution,  as  predicted  by  the  first  order  point  defect  model,  is  in 
good  agreement  with  experimental  data  (6).  Also,  while  three 
different  laws  (logarithmic  law,  inverse  logarithmic  law,  and  the 
parabolic  law)  are  necessary  to  explain  the  growth  characteristics  of 
the  passive  film  on  nickel  in  phosphate  solution  (6),  the  first  order 
point  model  is  able  to  account  for  the  growth  behavior  using  a  single 
growth  law  over  the  entire  range  of  conditions. 
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compared  with  those  inherent  in  the  empirical/semi-theoretical  logarithmic  and 
inverse  logarithmic  correlations.  The  purpose  of  this  comparison  is  to  assess 
how  the  various  integral  rate  equations  differ  from  one  another,  and  to 
determine  a  priori  whether  or  not  any  one  provides  a  more  acceptable  correl¬ 
ation  between  film  thickness  and  time  on  purely  mathematical  grounds. 

Logarithmic  Law 

This  semi -empirical  law  may  be  expressed  as  follows: 

L  =  A  +  B  •  lnt  (23) 

where  A  and  B  are  constants.  A  general  analysis  of  this  equation  reveals  the 
following  properties: 


(i) 

t  — ^  +  0 

lnt  — 

-  oo 

L 

->  - 

(if) 

t  - ^  +  oo 

lnt  — 

-f  CO 

L 

( i  i  i ) 

t  — >  1 

lnt  — 

+  0 

L 

^  A 

(iv) 

t  =  exp(-A/B) 

L 

=  0 

The  general  form  of  the  logarithmic  law  is  therefore  represented  by  the  graph 
shown  in  Figure  22.  The  most  notable  feature  of  this  plot  is  that  it 
resembles  that  predicted  by  Equation  (22)  at  times  greater  than  that  at  the 
inflection. 

Inverse  Logarithmic  Law 

The  inverse  logarithmic  law  is  usually  expressed  as  follows: 

1/L  =  C  -  Dint 

where  C  and  D  are  constants.  A  general  analysis  of  the  inverse  logarithmic 
equation  yields  the  following  characteristics: 
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Figure  21.  Schematic  film  thickness  versus  time,  as 
predicted  by  the  first  order  point  defect  model  [Equation 
(22)]. 


56 


€ x  I0“7, Field  Strength, Volt /m 


16 


P 

> 

► 


C0(L)=C0+(L)=CA.(0)=CA(0)sIOl4m"s 

■  €°«  IOa  Volt/m 

o  C0(0)=CD.(0)=CA.U.)=CA(L)=Ap=An= 


14 1—  IOMm-» 


L, Thickness ,  X 


Figure  24.  Calculated  electric  field  strength  vs  film 
thickness.  Note  that  the  points  designated  □  obey  the 
electro-neutrality  condition. 
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Figure  25.  Potential  function  A0(x)  vs  film  thickness. 

Note  that  the  points  designed  □  obey  the  electro-neutrality 
condition. 


Because  it  is  difficult  to  obtain  information  on  the  concentration  of  the 
different  charged  species  within  a  growing  oxide  film  and  at  the  two 
interfaces,  a  sensitivity  study  was  carried  out  to  determine  how  the 
calculated  space  charge  depends  upon  the  concentrations  of  various  species  in 
the  film.  The  various  cases  considered  are  as  follows: 

Anion-  and  Cation-Defect  Film 

In  this  case,  the  film  is  supposed  to  contain  all  six  charged  species 

V»  V(V*  v  v'»  v  /v  im»  e1,  h)  and  the  concentration  of  each  species  in 
u  Mx  Mu“u 

turn  at  each  boundary  is  varied  over  a  wide  range  (Figures  26  to  29). 
Anion-Defect  Film 

These  films  are  considered  to  not  contain  cation  or  ionized  cation 
defects.  The  sensitivity  of  the  calculated  space  charge  to  variations  in 
the  concentrations  of  anion  vacancies  at  the  metal /film  and  fi lm/ solution 
interfaces  was  investigated  (Figures  30  and  31). 

Cation-Defect  Film 

In  this  case,  the  anodic  film  is  considered  to  be  free  of  anion 
vacancies,  and  the  response  of  the  space  charge  to  variations  in  the 
concentrations  of  cation  vacancies  at  the  two  interfaces  was  determined 
( Figures  32  and  33) . 

The  electrical  neutrality  condition  can  be  expressed  by  the  following 
equation : 

-2Cd(0)  +  2Ca(L)  -  CD+{0)  +  Ca_(L)  -  An  +  AP  '  0 

If  the  equality  above  is  not  observed,  because  one  or  more  of  the 
concentrations  of  vacancies  becomes  greater  than  the  other,  the  strength  field 
through  the  film  does  not  remain  constant  or  close  to  constant,  but  changes 
abruptly.  Thus  space  charge  becomes  important  when  electrical  neutrality  of 
the  film  is  not  observed.  These  results  are  evidenced  from  Figures  26  to  33. 
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Figure  26.  Calculated  electric  field  strength  versus  film 
thickness.  Increasing  deviation  from  electrical  neutrality 
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Figure  27.  Calculated  electric  field  strength  versus  film 
thickness.  Note  that  in  this  case  the  anion  and 
ionized-anion  vacancy  concentrations  at  the  film/ solution 
interface  are  varied.  In  all  cases,  the  deviation  from 
neutrality  is  negligibly  small. 
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Figure  30.  Calculated  electric  field  strength  versus  film 
thickness.  Note  that  in  this  case  the  anion  and 
ionized-anion  vacancy  concentration  at  the  film/solution 
interface  are  varied  and  the  concentration  of  cation  vacancy 
at  interfaces  is  null.  In  all  cases  the  deviation  from 
electro-neutrality  is  small. 
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Figure  32.  Calculated  electric  field  strength  versus  film 
thickness.  Note  that  in  this  case  the  cation  and 
ionized-cation  vacancies  concentrations  at  the  metal/film 
interface  are  varied;  and  the  concentration  of  ion  vacancy 
is  null.  In  all  cases  the  deviation  from  electro-neutrality 
is  negligibly  small. 
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Figure  33.  Calculated  electric  field  strength  versus  film 
thickness.  Note  that  in  this  case  the  cation  and 
ioni zed-cation  vacancies  concentrations  at  the  film/solution 
interfaces  are  varied;  and  the  concentration  of  ion  vacancy 
is  null.  In  all  cases  the  deviation  from  the 
electro-neutrality  is  negligibly  small. 


[t  is  important  to  note  that  in  this  first  order  model  no  attempt  has  been 
to  relate  the  concentrations  of  various  species  using  equilibrium  or 
equilibrium  statistics.  For  example,  a  more  realistic  treatment  would 
ire  that  the  concentrations  of  electrons  and  holes  within  the  conduction 
valence  bands,  respectively,  be  related  to  the  concentrations  of  donors 
acceptors  by  the  Fermi-Dirac  distribution  function,  or  at  least  by 
zman  statistics  for  wide  band-gap  films. 

:ussion 


The  present  work  involves  the  calculation  of  space  charge  within  an  anodic 
n  based  on  a  model  that  involves  six  non-diffusing  species.  The 
culations  assume  a  steady-state  charge  distribution  within  a  film  of 
stant  thickness.  Accordingly,  a  growing  film  is  best  viewed  as  a  sequence 
steady-state  steps  in  which  the  film  thickness  increases  with  time.  The 
w  presented  here  is  that  of  one  such  step. 

This  model  is  different  from  that  developed  by  Hamilton  and  Miley  (4)  for 
iting  thickness  growth,  which  is  based  on  the  assumption  that  space  charge 
ses  the  bulk  concentration  of  the  rate  limiting  spaces  to  vary 
onentially  with  film  thickness.  In  the  present  case,  the  defect 
centrations  (donors  and  acceptors)  vary  according  to  Fick's  second  law  and 

electrons  and  holes  are  assumed  to  be  in  e .ectrostatic  equilibrium.  The 
centrations  of  these  latter  species  will  vary  exponentially  with  distance 
y  if  the  potential,  0,  is  a  linear  function  of  x.  Clearly,  this  occurs 
y  in  the  absence  of  space  charge.  Exact  solutions  for  concentrations  of 
rged  acceptors  and  donors  through  the  film  have  been  obtained  from  Fick's 
ond  law  under  steady-state  conditions,  and  a  numerical  approximation 
hnique  yields  solutions  for  the  dependence  of  the  space  charge  on  film 
ckness.  The  space  charge  is  calculated  subject  to  the  constraint  imposed 
the  Poisson-Boltzman  equation.  However,  the  film  as  a  whole  must  also  be 
ctrically  neutral,  although  to  our  knowledge  this  condition  is  seldom,  if 
r,  recognized  in  space  charge  computations.  In  this  particular  case,  the 
ctroneutrality  condition  has  been  applied  ex  situ,  although  in  later  work 
will  use  this  condition  to  fix  the  concentration  of  one  of  the  charge 
riers  at  a  boundary. 
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Sensitivity  studies  were  carried  out  in  order  to  determine  how  the 
ntial  through  the  film  can  be  affected  by  the  boundary  conditions  assumed 
he  model.  For  all  the  cases  studied,  the  space  charge  through 
film  was  found  to  be  imoortant  only  when  the  electrical  neutrality 
ition  is  not  satisfied.  However,  electrical  neutrality  is  a  physical 
irement  of  anodic  passive  films.  Accordingly,  for  steady  state 
itions,  it  is  apparent  that  space  charge  in  the  film  is  possibly  not 
rtant,  at  least  when  six  charged  species  are  included,  as  in  the  model 
idered  here.  It  is  important  to  note  that  the  boundary  conditions  (i.e. 
entrations  of  the  various  species  at  x=0  and  x=L)  were  varied  over 
emely  wide  ranges;  ranges  that  should  include  the  concentrations  of  the 
ies  at  the  boundaries  in  real  films.  It  is  believed  that  these 
ulations  help  to  clarify  the  controversy  which  exists  about  the  importance 
pace  charge  in  passive  films  (7).  Unfortunately,  comparable  experimental 
do  not  at  present  exist,  and  the  experimental  difficulties  involved  in 
rately  measuring  the  potential  drop  across  a  passive  film  have  not  been 
come.  However,  it  might  also  be  argued  that  the  steady-state  constraint 
d  simplistic,  and  that  the  lack  of  space  charge  computed  is  in  fact  an 
fact  of  the  model.  In  order  to  extend  these  calculations  to  include  the 
species  under  nonsteady-state  conditions  it  will  be  necessary  to  also 
Ider  the  generation  and  consumption  of  charged  species  within  the  film  as 
action  of  time  and  distance. 

IREAKDOWN  OF  PASSIVE  FILMS 

duction 


he  point  defect  model  that  was  previously  developed  (1,2)  provided  a 
stic  description  of  the  breakdown  of  passive  films  in  solutions 
ining  aggressive  anions  (e.g.  Cl”)  in  terms  of  the  accumalation  of  cation 
:ies  at  the  metal/film  interface  to  form  a  void  (cation  vacacancy 
isate)  of  sufficient  size  that  mechanical  collapse  ensues  to  nucleate  a 
The  enhanced  flux  of  cation  vacancies  across  the  film,  which  results  in 
accumulation  at  the  metal /film  interface,  is  considered  to  be  due  to  the 
:key  pair  reaction  and  the  adsorption  of  anions  (e.g.  Cl  )  from  solution 
inion  vacancies  (Vq*-)  at  the  film/solution  interface. 
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this  analysis,  the  point  defect  model  is  expanded  to  account  for  the 
of  minor  alloying  elements,  such  as  molybdenum,  on  the  anodic  breakdown 
live  films  on  metal  surfaces.  This  topic  is  important,  since  minor 
ig  elements  such  as  Mo  and  W  are  added  to  stainless  steels  to  improve 
• esistance  to  pitting  corrosion,  particularly  in  chloride-containing 
iments.  For  example,  the  addition  of  ->.2%  Mo  to  Type  304  SS  results  in 
L6SS,  which  is  frequently  used  for  seawater  service. 

Meal  Basis 

a  previous  study,  Chao  et  al .(1,2)  proposed  a  model  for  the  growth  and 
>wn  of  passive  films  which  yields  analytical  expressions  for  the 
il  breakdown  potential  and  the  incubation  time.  The  model  involves  (as 
t  approach)  transport  of  anions  (oxygen  anions)  and  cations  (metal  ions) 
ir  respective  vacancies.  Furthermore,  the  model  considers  that  anions 
?  from  the  film/sol  ution  (f/s)  interface  to  the  metal/film  (m/f) 
ace  and  subsequently  result  in  film  thickening  at  that  location.  On  the 
Hand,  cation  diffusion,  while  not  excluded,  results  only  in  dissolution 
t  in  the  growth  of  the  primary  passive  film.  This  model  has 
sfully  accounted  for  a  number  of  previously  unexplained  experimental 
ations  (1).  One  consequence  of  the  diffusion  of  metal  cations  from  the 

to  the  (f/s)  interface  is  that  metal  vancancies  (or  "metal  Holes")  are 
i  at  the  metal/film  interface.  These  metal  "holes"  will  tend  to 
je  into  the  bulk  at  the  metal  due  to  the  movement  of  cations  from  the 
into  the  film,  and  hence  to  "disappear".  However,  when  the  cation 
ion  rate  (i.e.,  the  metal  "hole"  production  rate)  is  higher  than  the 
f  "metal  hole"  submergence  into  the  bulk,  the  metal  holes  will  start 
up,  and  hence  will  form  a  void  at  the  metal/film  interface  (this  is  the 
s  of  pit  incubation).  When  the  void  grows  to  a  certain  critical  size, 
ssive  film  suffers  local  collapse  which  then  marks  the  end  of  the  pit 
tion  period.  Subsequently,  the  collapsed  site  dissolves  much  faster 
iy  other  place  on  the  surface  thereby  leading  to  pit  growth. 

:ording  to  this  argument,  the  criterion  for  pit  initiation  can  be 
sed  as  follows: 
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(jca  -  Om)(t  -  «)  >  5 

Where  is  the  cation  diffusion  rate  in  the  film,  J  is  the  rate  of 
submergence  of  the  metal  hole  into  the  bulk  metal,  t  is  the  time  required  for 
metal  holes  to  accumulate  to  a  critical  amount,  £,  (i.e.,  t  is  the  incubation 
time)  and  5  is  a  constant. 


From  the  first  order  point  defect  model  developed  earlier  in  this  report, 
and  considering  that  the  electric  field  strength  within  the  film  is  constant, 
Fick's  first  law  for  the  quasi-steady  state  may  be  solved  to  yield  the 
following  expression  for  Jca  (see  also  Refs.  1  and  2): 


Jca  '  xkDV 


Cv  (m/f)  exp(-XkL)  -  Cv  (f/s) 

*  Y  i  ¥  v  ' 

_jr _ mx 

exp(-XkL)  -  1 


where  X  is  the  charge  on  the  cation,  k  =  Fe/RT,  e  is  the  field  strength 
through  the  film,  F  is  Faraday's  constant,  R  is  the  gas  constant,  T  is 
temperature,  L  is  the  film  thickness,  is  the  diffusivity  of  the  cation 

V 

vacancy,  Cy  (m/f)  and  Cy  (f/s)  are  the  concentrations  of  cation  vacancies 
Mx  Mx 

at  the  metal/film  and  film/solution  interfaces,  respectively. 

As  previously  shown  (1,2)  it  is  possible  to  obtain  an  expression  for  Jca 
as  follows: 


J 


ca 


JV<X«> 


£x/2) 


where 


/n  \CW></2)] 

O'  =  xk£>„  —  •expt-AG"/RT) 

VW 
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/\\  /i(S(A-l)  '  FS-'!H  ‘  F«‘(f^)\ 

M  =  - „ - j 

n  =  Avogadro's  number 
n  =  molecular  volume  of  cations  in  the  film 
aG"  =  change  in  standard  Gibbs  energy  for  the  Schottky-pair  reaction 
aG®( a_d  =  change  in  the  standard  Gibbs  energy  for  the  reaction: 


VQ..  +  Cl"(aq. )  Cl q . 


ri0(f/s)  =  potential  drop  across  the  film/solution  interface  in  the  absence 

of  an  applied  voltage  and  for  pH  =  0 

a  =  activity  of  chloride  ions  in  solution 
Cl" 

a,B  =  constants 
Vapp  =  aPPlied  potential 

The  critical  pitting  potential,  Vc,  is  given  by  the  Equation  (36)  which 
leads  to 


1  -  nr(x/2)  pxn/  FaV^ 

Jm  ■ J  M  expNrr 


a-ci" 


(X/2) 


Therefore 


V 


c 


4.606RT 
X  Fa 


log 


2.303RT 

aF 


(37) 


and  the  incubation  time,  t.  is  a  function  of  breakdown  overpotential. 
Mathematically,  this  can  be  expressed  as 


^ind 


eyn/XFaA  V 
exp(  2rT 


'  -1 

1  +  T 
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It  is  well  known  that  minor  alloying  elements  can  have  a  profound  effect 
on  the  susceptibility  of  a  metal  or  alloy  to  localized  corrosion,  including 
pitting.  In  this  section  we  propose  a  solute/vacancy  interaction  model  ( S VIM ) 
to  account  for  this  phenomenon.  The  SVIM  is  develped  from  the  point  defect 
model  which  was  originally  proposed  to  account  for  the  growth  and  breakdown  of 
passive  films. 

The  basis  of  the  SVIM  is  that  minor  alloying  elements,  which  exist  in  the 
passive  film  in  a  higher  oxidation  state  than  the  host  cation,  can  interact 
electrostatically  with  cation  vacancies  as  described  by  the  following 
equilibrium: 


where 

K 

The  complex,  P,  and  the  solute,  A,  are  both  considered  to  be  fixed  in  space, 
and  only  the  cation  vacancy,  V,  is  assumed  to  move  through  the  film  from  the 
film/solution  interface  to  the  metal/film  interface. 


Assume  that,  in  the  absence  of  A,  the  concentration  of  cation  vacancies  at 


positon,  x,  is  ny.  The  concentration  at  the  point  x  after  complex  formation 


is  given  by  (nv  -  q[P]).  The  condi ton  for  breakdown  is  still  given  by 


(Jra  -  01.  (t  -  6)  >  £ 

ca  m  x=L 


in  which  the  flux  J  corresponds  to  that  for  the  uhcomplexed  vacancies.  When 

Ca 


no  complex  formation  takes  place 


\  /3n  \ 

7x=L=  V**/  Dknv(-6X) 


ca 


(38) 


but  when  complexes  exist,  the  flux  the  equation  must  be  modified  to  read 


H’  ■° 


a(n,  -  q[P]) 


ax 


+  Dk(-ex) (ny  -  q[P]) 


(39) 


where 


[P]-K„(nA-[P])(nv-q[P])c 


(40) 


and  where  nA  is  the  stoichiometric  concentration  of  solute  in  the  film 


Equation  (39)  can  be  writen  as  follows 


ca 


3n 


-0*  — —  +  0*k*(-ex)n. 


3x 


where  D*  and  k*  are  modified  to  reflect  the  presence  of  [pi .  These  functions 
must  be  calculated  for  specific  values  of  q,  i.e.,  for  various  complex 
configurations.  In  this  work,  only  q  =  1,  q  =  2,  and  q  =  3  are  considered. 
For  a  =  1  two  real  roots  exist  such  that 


where 


.  ,i/.!  \ 1/2 

»*/=-«*)/ 


D*  -  0  ?  1  1  -  nA"v)  1/2  (f  -  "A 


a  =  %  +  rnr +  nv 

Kn 


k*  =  Fe*/RT 


The  formulations  for  q  =  2  and  q  =  3  are  given  in  the  appendix.  Only  the 
real  roots  that  produce  positive  values  for  D*  and  k*  are  considered. 

In  order  to  quantify  the  effect  of  the  solute  on  the  pitting 
characteristics  of  an  alloy,  it  is  necessary  to  compute  the  changes  in  the 
breakdown  potential  and  induction  time  as  a  function  of  solute  (molybdenum) 
concentration.  The  appropriate  functions  are  therefore  given  by: 


/4.  606RT\  /  kD 

V  (XMJ  -  V  (XM  =0)  =( - log  - 

c  Mo  c  Mo  l  xFq  i  \k*D* 


tind  (XMo )  £  +  6  JV(x/2)[cn(x/2)  exp  J° 

AFaVc\" 

*ind  (XH0-°)  J‘  5  +  5  J'm"(x/2>CC1':,x/2)  exp 

This  last  equation  can  be  simplified  by  noting  that  for  6  — >  0 


tind(XMo) 

‘ind'Vr0’ 
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A  computational  algorithm  was  developed  for  an  Apple  lie  microcomputer  in 
order  to  obtain  the  difference  in  breakdown  potential  (Equation  41)  and  the 
ratio  of  incubation  time  (Equation  42)  for  different  concentrations  of 
molybdenum  when  q  was  assigned  values  of  1  to  3. 

Results  and  Discussion 


In  order  to  compute  the  effect  of  the  solute  on  the  pitting  potential  and 
induction  time,  it  is  necessary  to  specify  values  for  the  various  parameters 
in  Equations  (41)  and  (42).  The  parameters  of  interest  include  ny 
(concentration  of  vacancies),  nA  (concentration  of  solute),  and  «n 
(equilibrium  constant).  In  these  initial  calculations,  we  have  chosen  q=l  and 
the  parameters  Kn  and  n  are  selected  in  order  to  obtain  agreement  in  the 
first  experimental  point  in  the  plot  of  Vc  vs  molybdenum  content  for  18%  Cr 
ferritic  steel  given  by  Lizlov  and  Bond  (24)  (Figure  34).  This  graph  displays 
the  breakdown  potential  versus  concentration  of  molybdenum  for  pitting  of 
18Cr-%Mo  in  1  M  NaCl  at  25 °C.  Two  sets  of  parameters  resulted: 


20/cm3  ; 

1 

o 
•— » 

II 

c 

19/cn3  ; 

r- 

iH 

1 

o 

II 

Using  these  values,  it  was  possible  to  calculate  the  pitting  potential  (Fig. 
34)  and  ratio  of  incubation  times  (Fig.  35)  as  a  function  of  percent 
molybdenum  in  the  alloy.  The  generally-good  agreement  between  the  calculated 
and  experimental  breakdown  potentials  at  the  higher  molybdenum  contents  is 
most  encouraging,  even  though  the  "best"  values  for  q,  n  ,  and  Kn  are  so 
poorly-defined. 


A  sensitivity  study  was  carried  out  in  order  to  assess  the  importance  of 
each  of  the  parameters  in  determining  the  breakdown  potential  and  incubation 
time  according  to  the  SVIM.  Figure  36  shows  the  importance  that  the  value  of 
Kn  can  have  in  determing  the  pitting  potential  as  a  function  of  molybdenum 
content,  while  Figure  37  displays  the  influence  of  the  concentration  of 
vacancies  on  the  pitting  potential  and  incubation  time.  The  stoichiometry 
assumed  for  the  complex  also  exerts  a  strong  influence  over  the  passivity 
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Figure  35.  Effect  of  Molybdenum  on  induction  time  for 
pit-initiation  according  to  the  solute  vacancy  interaction 

model  (SVIM).  q-1,  nv=4xl020cm"3,  Kn=10'4cm3). 


VC(XMJ 


AV,  VxlO' 


%  Molybdenum 

Figure  37.  Influence  of  the  concentration  of  cations 
vacancies  on  pitting  potential  and  incubation  times  as  a 
function  of  molybdenum  content.  AV=V  (XM  1  -V  (XM  =0). 


Mo" 


breakdown  parameters,  as  shown  by  the  data  plotted  in  Figure  38.  Clearly,  in 
order  to  apply  this  model  to  the  prediction  of  the  effect  of  minor  alloying 
elements  on  the  pitting  characteristics  of  steels,  it  will  be  necessary  to 
generate  more  quantitative  data.  This  task  represents  a  major  challenge  in 
the  future  development  of  the  model. 
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APPENDIX 


CALCULATION  OF  BREAKDOWN  PARAMETERS  FOR  q  =  2  and  3 


The  general  equation  that  needs  to  be  solved  for  all  values  of  q  is  (Eq 


D  a£nv  " 


+  D  k(-ex)  (ny  -  q[P] j 


where  (Eq.  40) 


tp]-Kn  ["A  [",  -OflJ 


P  is  a  function  of  n.,  n  ,  K  andq.  n  ,  n.  and  K  are  constants  and  q  takes 

A  v  n  v  a  n 

the  values  of  1  to  3.  Equation  40  was  previously  solved  for  q  =  1;  here  we 
present  the  solutions  for  q=2  and  q=3. 


From  Equation  40  with  q=2  a  cubic  equation  in  P  is  obtained  as  follows 


3  2 

P  +  Pc 


(-V"vh  P  (vpA  *  "v2/4  -  •  Vy^4  ■  0 


The  solutions  for  this  case  are  then  developed  by  first  considering 


A.  =  -n.-n 
1  A  v 


A2  "  "vna  4  <"v  I74  -  1/<4Kn> 
2 

A,  =  -n.  n..  /4 


Q  =  3A2  -  A^ 


R  =  (  9AjA2  -  27A3  -  2Aj  )/54 


3  2 

D  =  QJ  +  R 


Thus,  the  roots  are  expressed  as 


P  -  S  +  T  -  Aj/3 


where 


S  =  (R  +  [Q3+R2]1/2  )  1/3 
T  =  (R  -  [Q3+R2]1/2  ) 1/3 


if  D<0  three  real  roots  are  obtained 


P  =  2  /^Q  •  COS  (8/3)  -  Aj/3 


P  =  2  v^Q  •  COS  (0/3+120)  -  Aj/3 


P  =  2  /=Q  •  COS  (0/3+240)  -  Aj/3 


where  8  =  arc  *  cos  (-R/(-Q3) ^2) 


P  =  Kn  (nA-P)  (ny-3P)' 


For  the  case  where  q  =  3,  Equation  35  can  be  written  as  follows: 


P4tP3  (-nA-nv)  *P2(pv2/3*nA-nv)tP(-PA-nv2/3-l/27Kn-nv3/27)  +Vr.#/27-0 


Defining 


8i =  vn. 


B2  ■  ",  «  *  nA% 


.  mm  \j 

•  .  -  "I 

1 

•  . 
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B3  ’  ‘"A", 73  -  1727  Kn  -  ",  727 

B4  =  ",V27 

and  noting  that  the  fourth  order  equation  in  P  can  be  reduced  to  a  third  order 
expression  by  the  following  substitutions 


A1  =  -  b2 


A2  “  Bl83  "  4B4 

A3  =  4B2B4  -  B32  -  Bj2  B4 


As  for  the  previous  case  (q=2)  it  is  possible  to  obtain  a  solution  for  P. 
Let 


M  =  S  +  T  =  Aj/3 


where 


R  =  (Bj2/4  -  B2+M)  1/2 


If  Rtf 


D  =  [3Bj2/4  -  R2  -  2B2  +  (ABjBg  -  8B3  -  Bj3  )  /4r]  1/2 
E  =  [38^/4  -  R2  -  2B2  -  (48^2  -  8B3  -  Bj3  )/4r]  1/2 


If  R=0 


D  =  ^3Bj 2/4  -  2B2  +  2(M2  -  4B4  j  1/2J 

1/5 

E  =  j3Bj2/4  -  2B2  -  2  (M2-4B4)  1/2J 


Thus,  the  four  roots  are  given  by 


P  =  -B  x/4  +  R/2  +  D/2 
P  =  Bj/4  +  R/2  -  D/2 
P  =  -Bj/4  -  R/2  +  E/2 
P  =  -Bj/4  -  R/2  -  E/2 

Once  values  are  obtained  for  [P],  for  q=2  and  q=3,  it  is  possible  to  write 
[P]  =  f  { n v ) 

but  nv  is  a  function  of  the  position,  x.  Accordingly, 

3[P]  9[P]  3ny 
3x  3n  ’  3x 

Let 


3[P] 


Therefore,  we  can  write 


3[P] 


3n. 


-  -  pp .  — _ 

3x  r  3x 


where  the  expression  PP  is  obtained  from  the  solution  for  [P]. 


Fick's  first  law  can  then  be  written  as  follows 


3n. 


\  v 

JcajL=-D  L1  -  PPJ  -  ^-*0k(-eX)(n 


-q[P]) 


where  PP  and  [P]  are  known  functions.  Letting 
D*  =  D  •  [1-PP] 
and 
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